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ABSTRACT 

We analyse the radial structure of self-gravitating spheres consisting of multiple inter- 
penetrating fluids, such as the X-ray emitting gas and the dark halo of a galaxy cluster. 
In these dipolytropic models, the adiabatic dark matter sits in equilibrium, while the 
gas develops a gradual, smooth, quasi-stationary cooling flow. Both affect and respond 
to the collective gravitational field. We find that all subsonic, radially continuous, 
steady solutions require a non-zero minimum central point mass. For mcgaparsec- 
sized halos with seven to ten effective degrees of freedom (F2), the minimum central 
mass is compatible with observations of supermassive black holes. Smaller gas mass 
influxes enable smaller central masses for wider ranges of F2. The halo comprises a 
sharp spike around the central mass, embedded within a core of nearly constant den- 
sity (at 10 1 — 10 2 5 kpc scales), with outskirts that attenuate and naturally truncate at 
finite radius (several Mpc). The gas density resembles a broken power law in radius, 
but the temperature dips and peaks within the dark core. A finite minimum tempera- 
ture occurs due to gravitational self-warming, without cold mass dropout nor needing 
regulatory heating. X-ray emission from the intracluster medium mimics a /3-model 
plus bright compact nucleus. Near-sonic points in the gas flow are bottlenecks to the 
allowed steady solutions; the outermost are at kpc scales. These sites may preferen- 
tially develop cold mass dropout during strong perturbations off equilibrium. Within 
the sonic point, the profile of gas specific entropy is flatter than s oc r 1 / 2 , but this is 
a shallow ramp and not an isentropic core. When F2 is large, the inner halo spike is 
only marginally Jeans stable in the central parsec, suggesting that a large non-linear 
disturbance could trigger local dark collapse onto the central object. 

Key words: accretion ■ cooling flows dark matter ■ ■ galaxies: clusters 
— hydrodynamics — X-rays: galaxies 



1 INTRODUCTION 

Galaxy clusters consist of baryonic and dark matter in 
the cosmic ratio (|Spergel et al.l 120071 ). Black holes and the 
stars in galaxies and in the intracluster light only consti- 
tute a small (10-15%) fraction of the baryons; and the in- 
tracluster X-ra y emitting hot g as comprises the ma jority 
of the baryons (|Lin et ajj 120031 iGonzalez et al.ll2007t ). Re- 
laxed clusters are found to contain a round core of approx- 
imately consta nt density, attenuating into fringes below de- 
tectio n limits (|Lea et al.l 1 19731 : ICavaliere fc Fusco-Femianol 
1 19761 ). There are also clusters with more centrally peaked 
core. They are thought to be systems with short radiative 
cooling time in comparison with the Hubble time. As radia- 
tive cooling causes the depletion of pressure support near 
the cluster centre, gas inevitably su bsides inwards from the 
cluster outskirts, i.e. cooling flows ICowie fe Binnevl [T977I : 
iFabian fc Nulsenll 19771 : iMathews fc Bregmanlll978l ). 



Cooling flows have been linked with accretion onto, 
and star formation in the dominant galaxy in t he cluster , 
and also the fuelling of th e ir galactic nuclei (e.g. ISilklll976l: 
Sarazin fc 0'Connell|[l983l IFabian et all 1984 iNulsen et all 



1984 ). Early models of cooling flow invoked a number of sim- 



ple assumptions. In some fluid formulations for the cluster 
structure, a static global gravitational potential was used, 
and there was no consideration of gas or halo self-gravity. 
Kinetic and ram-pressures were often not considered explic- 
itly, and this caused a cooling runaway near the cluster cen- 
tre, leading to a rapid deposition of a great amount of cold 
gas. Approximate deprojected cluster X-ray images indi- 
cated that t he gas inflow rates r h diminishes nearer th e clus- 
ter centre (|Stewart et al.lll984l : Frhomas et ail 1 19871 ). This 
was taken as evidence for widely distributed "mass dropout" 
- thermal instability spawning small, underpressured, invis- 
ibly cold clumps within a multiphase medium. Thermal con- 
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duction and magnetic fields were argued to be too weak to 
inhibit this instability and dropout. 

The early cooling flow models were challenged by 
various multi-wavelength observations (see r eviews by 
iDonahue fc Voitl 12004 iPeterson fc Fabianl [200^ ). First of 
all, radio and optical imaging have not shown the ex- 
pected accumulations of cold gas, nor the expected 
bursts of star formation. Moreover, X-ray spectroscopic 
and imaging deprojections of cluster profiles indicates a 
temperature floor typically a factor 3 or 4 below the 
peak temperature llKaastra et allfe oOl: Tam ura et al.ll200ll; 
Peterson et al]|200ll; Sakelliou et al.ll2002l ; Ijohnstone et al.l 
20021 ; IPeterson et al.1 I2003T ). In some systems the tem- 
perature even appears to increase at the smallest 

radii (e.g lO'Sullivan et al.1 120071 ). Spectral analyses sug- 
gest that the i ntracluster med i um (ICM) is like ly to 

be single-phase (|Bo hringcr ct al. 200 l| David et all l200ll ; 
iMolendi fc Pizzolatol l200ll ; iMatsushita et al.l I2002T ). These 
difficulties prompted the search of possible processes that 
could suppress the cooling flows, e.g. thermal conduction, 
or non-gravitational heating, such as th e power injected by 
activ e galaxies (AGN). (See review by IPeterson fc Fabianl 
120061 . and references therein.) There are still open ques- 
tions whether the heating processes can fine-tune to 
counteract the cooling stably, and whether the heating 
would distribut e appropriate l y acr o ss the relevant regions 
in the cluster |Fabian et al.l 1994 1_ Johnstone et al.l |2002| ; 



de Blok fc McGaughl 1 19971; IWeldrake et al.l 12001; Ide Blokl 



iBrighenti fc Mathewsil2003l ; IConrov fc Ostrikerll2008l ) 

Meanwhile, theories of halo structure have been over- 
turned several times. Once, it was assumed that cluster 
dark matter follows the distribution of galaxies, i n ap- 
proximately isothermal, flat-cored assemblages (e.g. iKind 
19661; iRood et all Il972); ICavaliere fc Fusco-Femianol Il97i 



Cowie fc BinnevHl977l ; (Fabian et al.lll98ll ). This view was 
naturally compatible with the classic signs (in the rotation 
profiles of disc galaxies) that galaxian dark matter is more 
shallowly and widely spread than the baryons. By the 1990s, 
cosmological N-body simulations were becoming fine enough 
to resolve cluster and galaxy halos, under the simplifying 
assumption that dark matter acts like a collisionless stellar- 
dynamical fluid (without any short- or long-range gauge 
fields of its own). Simulated halos develop sharp power- 



law central density cusps, (see e.g. r Dubiiiski fc Carlbere 
199l| ; iNavarro et "al]|l996l ; iMoore et all 19981 ; iDiemand et al 
2004 ; 



INavarro et all |2004| ; IDiemand et al.M2005l ). The 



distribution of coolin g, contracting gas may steepen the 
dark cusp further (e.g. Blumenthal et al . 1986: G nedin et al.l 
l2004ISelrwood fc McGaughll2005l ). 

In a circular way, cuspy profiles became an ansatz in 
the fitting cluster observations. Cuspy profiles have been 
assumed as templates in composite mass models fitted to 
gravitational lensing observations. It has been shown that 
a cuspy halo can hold a cored X-ray emitting gas distribu- 
tion qu alitatively similar to that of traditional cored halo 
models (|Makino et al.|[T998 | ). Gravitational lensing suggests 



flat cores in some clusters dTyson et al.lll998l ; Ig avazzi et al.l 
120031 ; ISand et all 12004 |200S| ). However, on galaxy scales, 
a considerable weight of evidence disfavours the existence 
dark cusps today (or implies that halos are less centrally 
concentrated than baryons). These lines of evidence in- 
clude ve locity fields of dwarf and low surface brightness 
galaxies (|Flores fc Primacklll994lMoorelll994l ; lBurkertll995l ; 



20051 ; ISimon et al.M2005l ; iKuzio de Narav et al.l 120061 ) ; and 
kinematics of dwa rf spheroidal galaxies ~l|Lokasl 120021 ; 
iGilmore et~alll2007l ). Current observational data for dwarf 
spheroidal galaxies ca nnot rule out cusp y profiles from 
the kinematics alone (|Walker et al.l 120071 ). though tidal 
tracers hint circumstantia ll y at gentle c ores in specific 
cases (|Klevna et al.1 120031 ; iGoerdt et all l200rj) . Lensing 
analyses of isolated elliptical galaxies suggest cuspy pro- 
files near the observed radii l|Read et al.l |2007| ) , while the 
kinematics of other cases imply fl a t cores or low dark 
densities dRomanowskv et al.l 120031 ; iDouglas et al.l 120071 : 
iForestell fc Gebhardt 20081 ) . Possible explanations of the 
cusp problem may involve subtle numerical systematics of 
N-body methodology, extra dark physics, or some forms of 
gaseous, stellar or AGN feedback. 

Feedback, if it is responsible for erasing cusps, must 
overturn a substantial fraction of a galaxy's baryonic mass, 
from the deepest zone of its potential, without leaving abnor- 
mal metallicities and stellar populations. The implementa- 
tion of feedback in numerical simulations suffers from severe 
challenges of resolution, and considerable arbitrariness or 
uncertainty in recipes for small-scale physics. Relevant fluid 
instabilities differ greatly between numerical schemes (e.g 
lAgertz et ai1l2007l ). Energy budgets of popular Lagrangian 
hydrodynamics methods are broken by en demic (but rarely 
mentioned) "wall heating" artefacts (e.g. lNohlll987l ). with 
unknowable consequences in simulated media where heat- 
ing, cooling or thermal instability are important. A defini- 
tive answer to the "feedback" question is far off; presently 
it is an almost unfalsifiable proposition. 

This paper aims to present a new formalism for the 
structures of relaxed galaxy clusters, and to probe the scope 
of its initial implications for cooling flows and dark matter 
(reserving empirical detail and observational fits for future 
refinements) . We reexamine the classic scenario of inflows in 
galaxy clusters with a more complete and consistent treat- 
ment of the gravitational interaction and energy exchanges 
in the gas and the dark matter components. Also a sensi- 
ble polytropic equation of state is used for the dark matter, 
that admits cuspless solutions for some systems and allows 
the multiple degrees of freedom in the dark matter. Note 
that a polytropic halo may arise if dark matter has strong 
self-interactions (SIDM), or if the system is formulated prop- 
erly in the framework of Tsallis' statistical the rmodynamics 
(cf. the Boltzmann statist ical thermodynamics. lTsallisill988l ; 
iPlastino fc Plastinolll993l ), or in collisionless systems with 
isotropic velocity distributions. We illustrate the properties 
and profiles of spherical, spatially continuous, stationary so- 
lutions relevant to cluster-sized systems. We quantify certain 
signature radii of these solutions, for the benefit of compar- 
ison with simpler models in the literature, and to inform 
future observational tests. Our solutions indicates that it is 
inevitable that point-like central masses would emerge in 
relaxed clusters, groups or pressure-supported galaxies. For 
some appropriate regimes of the halo micro-physics, the pre- 
dicted minimum central mass is consistent with those of su- 
permassive black holes in giant galaxies. We discuss implica- 
tions for the rapid origin of supermassive black holes, mono- 
lithic condensation of early stellar populations in galaxies, 
and the problem of the central structures of dark halos. 

The paper is organised as follows: §2 gives the general 
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formulation of the multi-component self-graviating systems, 
and the construction of the gas and dark-matter dipoly- 
tropes. Next, §3 describes the valid solutions in the cluster 
parameter-space, the properties of the solutions, their com- 
parisons with current observations, and predictions to be 
tested by future observations. In §4 we discuss our results in 
the context of galaxy and cluster evolution. We conclude in 
§5. The appendices show the derivation of our model's nat- 
ural units, the normalisation and the rescaling of the model, 
the interpretation of the effective degrees of freedom for the 
dark matter and comparisons between our model and other 
standard spherical models for clusters. 



2 MULTI-COMPONENT SELF-GRAVITATING 
SYSTEMS 

2.1 multi-fluid formulation 

In our model the system has multiple components. It is self- 
gravitating, and the distinct components interact among 
each other through their shared gravitational potential. 
Each component, i, has its equation of state, which takes 
the form 



Pi = piOi = SiPt 



(1) 



with partial pressure pt, density pi, isotropic velocity disper- 
sion ai (which corresponds to an isothermal sound speed in 
the fluid description), specific entropy s* and adiabatic in- 
dex 7;. The adiabatic index is related to an effective number 
of degrees of freedom Ft via 



7> = l + ? r 



(2) 



-pi + V ■ piVi = , 
-/XVj + V • piViVi + Vpia't = pdi , 
-e» + V ■ (e» + piot)vi = pw -fi + d , 



(3) 
(4) 
(5) 



The mass, momentum and energy conservation equations 
read 

d_ 

of 
d_ 
at' 
d 

di L 

where the energy density is 

_ Pi a i i 1 2 S t p]' 1 2 fe .s 

e ^ — +2 P ^ = — + 2 P ^ ■ (6) 

Equivalently, the energy conservation equation ([5]) may be 
expressed in term of the entropy: 

+ Vi • Vsi = (ji - l)dpp' ■ ( 7 ) 

The variable d is a volumetric power, which specifies the 
energy gains and losses. For a component net loss, say ra- 
diative cooling, d < 0. We assume that d is determined 
by the local thermodynamic and dynamic variables. This 
assumption is justified in the cluster environments as the 
radiative processes are optically thin. 

The gravitational potential, $, satisfies the Poisson 
equation 



V 2 $ = AnGj^Pi 



(8) 



The gravitational force is determined from the gravitational 
potential f = — V3> for all i, and the force field is the same 
for all the mass components. 



2.2 steady spherical solutions 

In this study we consider only stationary spherically sym- 
metric systems. Time dependent analysis will be discussed 
in Saxton et al. (in preparation), and systems with more 
complicated geometries will be presented elsewhere. In a sta- 
tionary spherically symmetric system, the density, velocity, 
temperature and gravitational field are functions of radial 
co-ordinate r only. The mass continuity equation ((3} requires 
that density and velocity are related by 



rhi = 4nr piVi 



(9) 



For mass inflow rhi < 0; and for mass outflow rhi > 0. A 
system with a hydrostatic equilibrium has rhi = and Vi = 
everywhere. The entropy equation (0 becomes 

dSi f * \ Li 

rf7 = (7i 



ViPi 



(10) 







and the conservation equations (J3j — dSJ) now read 
Pi 

piVi pi 
Pi v i PiVi 

where the three source/sink terms are 
„ 2piVi 

■Zli = , 



dr 



pi 




Z\i 


Vi 




Z^i 






Zzi 



(11) 



(12) 

(13) 
(14) 



Zn = Pif , 
Z-a = piVif + d 

The component mass m; interior to radius r is given by 

drrii 



dr 



■ 4nr pi , 



and is related to the force by 



It follows that 

df 2/ ^ 

— = 4nG > pi 

dr r t— 1 



(15) 



(16) 



(17) 



For a system with a central point-mass (e.g. a supermassive 
central black hole of a cD galaxy in a cluster) / rises asymp- 
totically near the origin. For systems without a point mass 
(323, / = at the origin. 
Inversion of (|ll|l gives 



dpi 
dr 

dvi 
dr 

da[ 
dr 



1 

A~ 

1 

A~ 



+ Pif- 



(7< - 



27i CT; Vi 



«i/ + (7i-l) — 

Pi. 



Ai 



o„.2 2 2 
^Vi <T; 2 , <7i 

h ff» / + ■ 



v 



PiVi 
,2 „,2 



(18) 
(19) 
(20) 



with the sonic factor, Ai = "ficr? — vf, where Ai > corre- 
sponds to a subsonic flow. If the matrix in (|ll|l is singular, 
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one of the hydrodynamic variables must be eliminated alge- 
braically to reduce the number of equations. For stationary, 
spherically symmetric systems, two of (pt,Vj,crf ,pi, si) suf- 
fice to describe the stationary solution. 

We will distinguish the mass components in the model 
for zero and nonzero rhi, as the solution to above equations 
depends on whether or not the mass component is in a bulk 
inflow. In either case the mass conservation equation Q will 
be used to eliminate one of the dynamical variables of each 
component. 

2.3 stagnant component, rhi = 

Usually, a static structure with zero inflow (rhi = 0) is for- 
bidden (see Equations 1181 and I20p . if there is a net energy 
loss (d 7^ 0) or if the system is not isothermal (7, 7^ 1). An 
isothermal condition can be established if heat transport is 
more rapid than other radiative and dynamical processes 
in the system. This is not easily satisfied for the radiative 
gas components in a cluster. However, the situation is dif- 
ferent for the dark matter component as dark matter nei- 
ther radiates nor absorbs light. Self-interacting dark matter 
could behave like a fluid. Thus, such a dark halo may have 
a hydrostatic profile satisfying rhi — 0, d = 0, Vi = 
and dvi/dr — everywhere. Moreover, its structure is com- 
pletely specified by the density {pi) or temperature (velocity 
dispersion) (of) profiles. 

The gradients of the density (|18|l and temperature (|20[) 
are given by 

dpi _ Pi I 
dr ji of 

and 

d<Jj _ Ji - 
dr 7, 



(21) 



(22) 



respectively. The latter implies that a 2 = — (7* — 1)($ — 
$fl)/7i- In the cluster setting, corresponds to the dark 
halo surface potential. 

The specific entropy is uniform in the cluster. The ef- 
fectively polytropic equation of state is appropriate for dark 
matter if it is self-interacting, governed by Tsallis thermo- 
statistics, or has a power-law phase-space density (i.e. col- 
lisionless). In the latter case, the momentum equation is 
identifiable as the Jeans equation in the isotropic limit. (See 
appendix [C]) We note that if only one self-gravitating fluid 
component is present and if there is no central point-mass, 



the object is essentially a Lane-Em den sphere l|Lanej fl870l : 
lEmderj|l907l ; IChandrasekhaill 19391 ) . 



2.4 flowing components, rhi 7^ 

Fluids with m, / exhibit a central density cusp, p; — > 
00 as r — > 0. In order to keep the equations numerically 
tractable, we use ((9| to eliminate pi. The equation for the 
density gradient (|18[1 is then redundant, leaving two relevant 
gradient equations: 



dvi 
dr 

da[ 
dr 



f + 



(7. ~ 1) 



m , 



in - 1K 2 



/ + — + 

r 



(23) 



(24) 



and the mass profile equation: 

drrii rh 
dr Vi 



(25) 



Elimination of pi in terms of m and r, particularly from 
the cooling function d, clarifies the asymptotic behaviour 
of the differential equations, especially in the inner regions. 
The flow velocity Vi may take any value at the origin. Mod- 
els in which Vi 7^ at the origin describe accretion onto, or 
winds emerging from, a central object, presumably compact. 
Note that inflowing constituents have some cuspy behaviour 
near the origin (either p; — > 00, Vi — > or Vi — > 00) regard- 
less of the occurrence of radiative cooling. This causes some 
trouble in numerical integration. Changes of variables, in- 
cluding a switch of the independent variable, may, however, 
circumvent the problems. 



2.5 composite system: radiative gas embedded in 
self-interacting dark matter 

We consider a model cluster, which is a composite system 
consisting of a cooling gas component with rhi < and 
£1 < and a self-interacting dark matter component with 
rh2 = and £2 = 0. We treat both of them as separate 
fluid components. We omit the stars in galaxies, as they 
are a minor fraction of the cluster baryons. The dominant 
radiative loss of the gas component is optically thin ther- 
mal bremsstrahlung radiation. We omit line cooling, dust, 
conduction and Compton effects. Thus, the radiative loss is 
specified by a cooling function: 



£1 



-Bp\Ot = —Bp\ 



(3+7l)/2 1/2 



(26) 



wher e the normalisation B dep ends on the gas composi- 
tion l|Rvbicki fc Lightmanl 1 19791 ). Cooling domination im- 
plies Ci < everywhere, and the specific entropy increases 
monotonically with r. This also ensures buoyant stability. 

The inner boundary is a gas density cusp, pi — > 00, re- 
gardless of whether radiative cooling or compressional heat- 
ing dominates the gas inflow (accretion). There are two 
kinds of cusps: the cold cusp (cti — > 0) and the hot cusp 
(<Ti — > 00). In cold cusps, the gas pressure pi is finite every- 
where. The specific entropy vanishes, si — ► 0, at the cusp. 
For the hot cusps, the entropy decreases smoothly towards 
the origin. The flow reaches the origin exactly, r* = 0, where 
the speed and pressure reach infinity. For fluids with an adi- 
abatic index of 7; = 1 + 2/ Fi, the asymptotes are p oc r~ Fi ^ 2 , 
a 2 oc r -1 , and v oc r ( Fi ~ 4 ^ 2 . There may be a positive-mass 
compact object at the origin, m*. The flow at a hot cusp is 
a self-gravita ting generalisation of subsonic Bondi accretion 
(|Bondilll952h . 

A fully general model of multi-fluid self-gravitating ob- 
jects permits supersonic inflows, existence of sonic points 
and formation of shocks. However, we focus on quiescent 
systems in this study. We consider solutions in which any 
inflow, if present, is subsonic everywhere in the cluster. 



2.6 central asymptotic behaviours 

The asymptotic power-law behaviours of the variables in the 
central region allows us to introduce a new set of variables: 



Pa 



(27) 
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Pvi 



Pi r 



Fi/2 
„(4-fi)/2 



(28) 
(29) 



where i = 1 for the gas and i — 2 for the dark matter. These 
variables are finite at r = 0. The corresponding equation of 
state is given by 



PtJi — ^i-Ppi 



(30) 



We define a logarithmic radial coordinate I = lnr. In terms 
of the new variables, the gradient equations are now 



d/3 Vi 
dl 



dl 



d/3 Pl 
dl 



Pvi 



i-Fi 



-Gm 



2fiPc 

2 
F 



Plt c 



,(31) 



= P° i + 



7,-1 



= P Pi 



F 
2 



7;/Ml-X 2 ) 



2 7l /?a,A4 2 



-Gm - (1 - 7l M 2 )/3ir c 



2 7i /3 CTl >t^ 



(32) 



—Gm + —Plt c 



where m = mi + m 2 . The entropy equation is 



-rr =-8i< (7i - 1 )= 



The cooling function 



3l 



Pv 1 



and the radial index of cooling term 
_ 7 - 2Fi 



,(33) 



(34) 



(35) 



(36) 



The Mach number M 2 = fi 2 vl r Fl 3 / 7 i/3 CT1 , and its profile is 
given by the equation: 



dl 



1-M 2 



-(1+7^)^^ 
whose solution to the profile equation requires that 



IimAr 



oo if Fi < 3, 

Ml > if Fi = 3, 
if Fi > 3 



(37) 



(38) 



Note that the dark matter component does not have 
radiative cooling. This implies (3 V2 = 0, d/3 V2 /dl = and 
dsifdl — throughout the system, and the dark- matter 
structure is determined by 

2Gm 



dl 




d/3 P2 
dl 


= ( F i 



72-1 

72 



Gm = f3 a 



F 2 + 2 



(39) 



a _ Fi R (F 2 -2)/F 2 d(3 a2 

7 2/3 CT2 y^ 2_ 2s/ p - di ■ (W) 



The mass and moment of inertia are given by 

drm a (e-FA/2 



and 

dh 
dl 



8 Ti 



Pp 



(10-^)72 



(41) 



(42) 



respectively. The mass profile would have a steep, cuspy gra- 
dient near the origin for Fi > 6. The moment of inertia shows 
a central cusp when Fi > 10, but this Fi corresponds to sys- 
tems with infinite mass and radius, which are unphysical 
and are irrelevant to astrophysical galaxy clusters. 



2.7 numerical calculations 

Some numerical difficulties could arise in solving the struc- 
ture equations given in the above section when Fi > 6. To 
overcome these we consider another coordinate 



„(6-F 2 )/2 



6- F 2 

instead of I and the transformation 

= r .(F 2 -6)/2 = 2 

da (6 — F2)a 



(43) 



(44) 



for the derivatives. 

We set the boundary conditions at the surface of the 
dark-matter component, R, and the integration will pro- 
ceeds inwards to the cluster centre. At the outer bound- 
ary (r = R), we specify the total mass (m(R) = mi(R) + 
m2(R)), the matter inflow rate (m), the specific entropy 
of dark matter (a constant s 2 > 0), the gas temperature 
(Tr = cr'i(R)) and Mach number (Mr = M(R)). The den- 
sity and temperature of the dark-matter are zero at r — R. 
The entropy s 2 is non-zero, and it defines the dark-matter 

density and temperature gradients. 

An adaptive-step Runge-Kutta scheme (|Press et al.l 

Il992|) is used in the integration. We first integrate a small 
step radially inwards, using a 2 as the independent vari- 
able, to avoid numerical troubles that could be caused by 
the steep gradients at the cluster boundary surface. We 
then proceed with the main integration, using the vari- 
able I or a, approaching a reference radius chosen to be 
r = lO -15 ^. (Here U x = B/G is the natural unit of dis- 
tance; appendix [X]) We examine how the variables behave 
near this radius. If the gradients become too steep, we would 
consider alternative variables for the integration. There are 
two types of breakdown that may necessitate a switch. We 
name them as "cold catastrophe" and "supersonic catastro- 
phe", and will discuss each of them in more detail. 

The "cold catastrophe" arises when the cooling is too 
efficient, causing the temperature to plummet steeply. We 
make use of the local monotonicity of si to define a new vari- 
able z = s\^ 2 for the integration. Although there are a steep 
radial gradients for the variables f3, i.e. the corresponding 
\df3/dl\ diverge, the derivatives dfijdz are still well-behaved 
and finite. Thus, it allows a smooth integration towards the 
centre, where z — ► monotonically. Note that if the cooling 
catastrophe occurs at a non-zero radius - forming a zero- 
temperature shell - we may discard it as it is not a viable 
steady solution. The cold shell lacks pressure support, and 
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material at large radius would fall inwards until a more sta- 
ble configuration emerges. Such a system would have a va- 
riety of interesting dynamic behaviours, and we will discuss 
it and related systems in a separate paper. 



The "supersonic catastrophe" arises when the gas Mach 
number increases towards unity at a certain radius. The 
(1 — M 2 )^ 1 factor will diverge and create numerical prob- 
lems. As a resolution, we switch to M 2 as the independent 
variable when the integration proceeds and approaches the 
sonic horizon (M 2 — > 1). Two situations would occur. In 
the first one, the solution has a discontinuity, with causal 
disconnection between the inner and outer regions. This is 
the shock solution, which is interesting but does not corre- 
spond to steady galaxy clusters, the prime interest in this 
work. The second one correspond to a smooth tr ansonic 
flow, in which the inflow would pass a sonic point (|Bondil 
1952), beyond which the accretion becomes supersonic. The 
steady transonic solutions are valid, but they tend to give 
lower gas densities than the solutions with subsonic inflows 
throughout the entire cluster. This also implies a greater 
residual central mass m* . In this paper we seek to minimise 
m„ within the set of truly steady solutions, and we prefer 
the wholly subsonic solutions. 



We consider various trial (Mr,S2) at R in the inte- 
gration. In each trial we record the radius, r* where inte- 
gration stops, and the central, interior mass, m„. For fixed 
{Fi, F2, R, m(R), rh, Tr), we map the (A4r,S2) plane and 
divide it into zones according to their physical and numer- 
ical characteristics (see Figure [TJ. Qualitatively, we have 
four principal zones. Three of them are either unphysi- 
cal or irrelevant to astrophysical galaxy clusters: (i) "too 
cold" — afflicted by a cooling catastrophe at a certain ra- 
dius; (ii) "too fast" — containing a supersonic disconti- 
nuity; (iii) "levity" — with insufficient pressure support, 
implying a negative central gravitating mass in compensa- 
tion. The acceptable, physical solutions lie in the wedge- 
shaped region between the "too cold" and "too fast" zones. 
These solutions correspond to steady structure and subsonic 
flow throughout the cluster. The {Mr, S2) values where the 
"too cold" and "too fast" boundaries intersect depend on 
(Fi, F2,R, m(R),rh, Tr). The tip of this wedge region of the 
acceptable solutions does not generally touch the contour 
where m» = 0, i.e. the boundary of the "levity" zone. The 
central mass m* of a physical cluster must be positive and 
exceed some certain value. All steady, self-gravitating, spher- 
ical, cooling multi-component clusters would require a cen- 
tral point mass. Strictly speaking, a spherical, cooling clus- 
ter with gas and dark matter composites without a central 
mass condensation is never steady. It will eventually evolve 
into another configuration on a dynamical timescale. Read- 
justment may start with the growth of a central mass, with 
waves of disturbances propag ating outwards like t he "s wal- 
lowing waves" as described in iMathews fc Bakerl (| 1 9 T lh . In 
the later sections, we will present two of the mechanisms 
that lead to the formation a central condensation in multi- 
component galaxy clusters. 
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Figure 1. Broad and detailed maps of the key parameter do- 
mains of models with standard mass m(R) = A0U m ~ 3.57 X 
10 14 7710, radius R = 4U X ss 0.983 Mpc, inflow m = 10m Q yr" 1 , 
gas surface temperature 1 keV and degrees of freedom Fi = i*2 = 
3. The axes are the surface gas Mach number and the dark specific 
entropy. In a wedge-shaped domain ("deep", top panel) steady 
models reach from the dark surface to the origin. To the right 
("too fast") the profile suffers a supersonic break at an interme- 
diate radius. To the left ("too cold") a cooling catastrophe occurs. 
Models in the shaded region ("levity") require a negative central 
mass. The "fast" and "cold" borders intersect at a point (circled) 
above the m* = contour. Thus m„ has a positive minimum for 
truly steady models. 



3 STEADY INFLOW SOLUTIONS 

3.1 size, mass and compositional families 

We now compare model clusters with different inflow rates 
(rh), gas surface temperatures (Tr) and dark- matter degrees 
of freedom (F2) for a given total mass m(R). We choose 
to fix m(R) = 40U m ~ 3.56 x 10 14 m©, unless other- 
wise specified, although the masses in the cluster solution 
are rescalable (see Appendix [Bj. We fix Fi —3 for the 
gas. The values of rh span the range 1 m© yr _1 ^ rh ^ 
1000 m© yr -1 , inferred from X-ray imag ing and spectral ob- 



servations of cooling core clusters (e.g. iFabian et al. 
Nulsen et all 1 19821; Istewart et all Il984l ; IWhite et all 



1981 



1994 



Edge et al"1ll994l ). For given (Fi, F2, R, m, rh, Tr), we min- 
imise m* over the (Mr,S2) plane. The contours of m» and 
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gas fraction (1/T = mi/m) are plotted in the (F2,R) plane 
(Figure [2]), with other parameters held constant. 

The radius of a minimal-m, cluster increases with F2 
along a T-track, All else being equal, gas-richer tracks have 
larger cluster radii. Each family of solutions shows a simi- 
lar variation of R with F2 when following a particular T- 
track: for cosmic composition, we find that R(F2 = 9) fa 
1.229i?(i<2 = 2). Also, m* decreases as F2 increases in a T- 
track. However, for fixed F2, m* varies with R. As shown 
in the upper panel of Figure [2] a peak m» appears for 
R fa 1.25 Mpc, in the models with F2 ~ 8, rh = 1 ttiq yr" 1 
and T R = 1 keV. 

The m* contours behave qualitatively differently in sev- 
eral distinct regions of the (F2,R) plane: 

(i) For F2 > 7, the values of m, drop steeply with increas- 
ing F2. In the case of rh = 1 yr _ , m* drops by a factor 
~ 0.1 for each increment of 1 in F2. The drop is steeper for 
larger rh. The cooling and sonic constraints permit a smaller 
central mass if the central density profile is steep. This oc- 
curs most easily for halos with more degrees of freedom. 

(ii) In another regime, with small radius R, both the gas 
fraction and the minimal m* drop steeply with decreasing 
R, regardless of F2. The poverty of gas loosens the "cold" 
and "fast" constraints, enabling smaller m„. As R shrinks, 
the solution approaches that of the Lane-Emden ideal poly- 
trope (which lacks a central mass) or else it becomes a point 
mass lacking both halo and gaseous envelope. 

(iii) The rest of the (F2,R) plane is a relatively feature- 
less plateau (top-left region of Figure [2} : m* increases only 
slightly even when there is a large increase in R. The T- 
tracks, however, vary smoothly across the borders from 
plateau to the low-m* slopes. This insensitivity of T occurs 
because the determination of bulk composition is global, 
whereas m* is governed by local gas constraints acting in 
local bottlenecks of the inflow at small radii. 

The attainable range of m, values across the (F2 , R) 
map depends on the gas inflow rate, rh. Smaller rh reduces 
variation in m„, with lower values on the plateau. The three 
panels of Figure [5] compare families of solutions that differ 
only in rh. The m, contours are completely different for 
different choices of rh. For rh = 1 m@ yr _1 , the cosmic- 
T track crosses ~ 4 orders of magnitude in m„. For rh = 
10 mo yr _1 , the equivalent track crosses almost ~ 4.5 orders 
of magnitude. For rh — 100 m© yr _1 , the track crosses ~ 5 
orders of magnitude. Note that increasing rh also shifts each 
T-track to smaller radii, i.e. for a given composition, clusters 
with strong inflows tend to be more compact. 

Increasing the gas surface temperature Tr (with every- 
thing else fixed) shifts the T-tracks to smaller radii. However 
varying Tr has negligible effect on the m* contours. Thus 
the m* — m*(F2) profile of a given track shifts to slightly 
lower masses, but this is only due to migration of the track 
across fixed m„-contours. (Compare the lower panel of Fig- 
ure [5] with Figure [3] Thus, all else being equal, a hotter 
cluster is a smaller cluster but with a similar central mass. 
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Figure 2. For a fixed cluster mass (m = 40), surface tempera- 
ture (Tr = 1 keV) and inflow rate, we vary the outer radius R and 
dark degrees of freedom F%, Black/solid contours map the mini- 
mal values of log 10 (m* /m,Q), the central point mass. Red/dashed 
contours show the gas fraction, 1/T relative to the cosmic baryon 
fraction (i, i ~ , 1, | . . .)/Too. This sequence of panels shows the 
effect of varying the inflow rate, with rh = 1, 10 and 100 ttiq yr~ 1 
from top to bottom respectively. For large rh, the T-tracks occur 
at smaller radii, and m* varies more widely along each track. 



3.2 external atmosphere 

As rh is a constant at all radii, all stationary solutions have 
some tenuous gas extending indefinitely far beyond the dark 
halo. This atmosphere must lose its identity in the Hubble 
flow at some point. In this work we take this cosmic atmo- 
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T R =5 keV, rh = 1 00 m yr 1 F 2 =8, rh =1 m yr 




Figure 3. Minimal-m* and 1/T map corresponding to the bot- 
tom panel (100 ttiq yr^ 1 ) of Figure |21 but for a warmer (3keV) 
cluster surface. The m, contours are almost unchanged. The 
tracks for given 1/T occur at smaller radii. 
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Figure 4. Minimal-m* and T-track map as in Figure[2] but with 
hotter gas and more inflow (Tr = 6 keV, m = 1000 rag yr ). 
This map includes our most compact solutions of a given T. 



sphere or accretion flow for granted. We will not consider 
its distribution in detail, but briefly discuss the qualitative 
implications of two scenarios. 

We may suppose that infinite atmosphere is a contin- 
uation of the cooling, polytropic gas inflow but without a 
dark component. The density attenuates with distance, and 
bremsstrahlung cooling becomes negligible. The effects of 
local self-gravity dwindle. If no sonic horizon appears, then 
the asymptotic atmospheric structure follows some power- 
law decline. 

Alternatively, we may choose to in terpret the halo sur - 
face gas as a cosmic accretion shock (e.g. lBagchi et all2006l ). 
Its temperature depends on the cluster mass and radius. In 
practice we select m(R) and Tr. Then shock conditions con- 
strain the plausible range of radii, 



Fi 



Gm 



(F a + l) 2 T a 



Fi Gm 
Fi + 2 ~Tr 



(45) 



with the lower and upper limits corresponding to strong and 



Figure 5. Minimal-m* and 1/T map for F2 = 8, m = 
10 mg yr - 1 and variable surface temperature Tr. The T-tracks 
shift, but the minimal m* hardly changes. 



weak shock extremes respectively. More exact constraints 
emerge if we consider the post-shock Mach number and pre- 
shock cooling. Not all families of models enable the cosmic 
T-track to satisfy (|45fl . For rh = 1 m© yr" 1 and Tr — 1 keV, 
the radii of (to*, T)-optimal solutions are too large, however 
a cooler family of models with Tr — 0.4 keV is satisfactory. 
By adjusting m and Tr in the opposite direction, we obtain 
hotter and more compact cluster models that also suit a 
shock interpretation (e.g. Figures [3] and [4} . 

For very high Tr the atmosphere would extend in- 
finitely, with pressure dropping to some asymptotic value, 
and temperature rising as a power-law. This restricts 
the upper limit of Tr for physical solutions. Such phe- 
nomena also occur in t he adiabatic, hydrostatic clusters 
(|Gull fc Northoverlll975h . where below some threshold (ef- 
fectivel y a minimum Tr), a gas i nflow truncates at finite 
radius (IMathews fc Bregmadll978h . 

In our solutions, flow continuity (constant to) means 
that the sub-critical gas atmospheres can break at an ex- 
ternal sonic point {M 2 = 1) with non-zero density. It is 
unclear what external conditions should match onto such a 
supersonic break. Reducing Tr further causes the gas pro- 
file to break somewhere inside the halo, r < R. These are 
not numerically feasible, searchable solutions, and we avoid 
them. We have coarsely scanned the parameters (Tr, R) for 
fixed (F2, m) and found little qualitative variation in the in- 
ner profiles or minimal-m* values, aside from the T-tracks 
shifting (see Figure [5j. 



3.3 radial structure of particular clusters 

3. 3. 1 general properties and density profiles 

We now examine the internal structures of specific clusters 
in detail. Here we discuss only the minimal-m* models where 
the overall composition is cosmic, 1/T w 0.163. Table[T]lists 
the parameters and some global properties of these clus- 
ter solutions. In each model we chose F2 for the halo and 
(rh, Tr) for the gas, then tuned the cluser radius R to obtain 
cosmic composition. We tabulate signature radii of the mod- 
els defined in Appendix IDll Ri 1 , Ri 2 , Ri are effective core or 
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lever radii, weighted by inertial moments of gas, the halo and 
both combined; R w characterises the concentration of gravi- 
tational potential energy; Ri, R2, Rs,Ri are radii where the 
total density has a radial logarithmic slope of —1,-2,-3 
and —4; R is the outermost peak of the rotation curve. 
Appendices [D2HDTT1 compare these signature radii to other 
spherical models in the literature (see Table [DTI for data). 

The gas + halo models differ from gasless models (Ap- 
pendix ID2f) in several key respects. For F2 > 7 the gassy 
clusters are gravitationally more compact: R m /R is smaller 
than for corresponding gasless, non-singular polytropes. The 
central mass and dense cusp deepen the potential well sig- 
nificantly (especially when F2 is large). 

The concentration of gas mass (Ri ± /R) is rather insen- 
sitive to F2. In all of our minimal-m* models, the gas is 
less centrally concentrated (0.74 < Ri t /R < 0.81) than in 
a simple F = 3 polytrope (Ri/R ~ 0.715). The presence 
of gas affects the halo concentration (Ri 2 /R), depending on 
i*2. For F2 = 2,3, the dark mass is slightly more concen- 
trated (smaller Ri 2 /R than for gasless spheres). For larger 
F2 , gas makes the halo less concentrated (larger Ri 2 / R than 
for polytropes). The combined mass distribution has an ef- 
fectively intermediate concentration: either Ri 1 < Ri < Rj 2 
or Ri 2 < Ri < Ri 1 . 

The rotation curve peaks farther out than in gasless 
halos (Lane-Emden spheres) of the same F2, and moreso for 
large F2. Specifically, R enlarges by ~1% for F2 = 2 but by 
~20% for F2 = 9 cases. In all cases we find that R a > R3. 
Importantly, this means that the rotation curve doesn't peak 
until outside a radius where the combined density slope is 
steeper than —3. For optimised gassy models with F2 = 2, 3 
we find R > R4; however we find R3 < R < R4 for F2 J? 
5 generally (and for F2 ^ 4 for the compact models). In 
contrast, the gasless halos have R2 < Ro < Ri for F2 = 
2, 3, 4, 5. For real, relaxed galaxy clusters, the comparison of 
measured R3 and R4 (e.g. from gravitational lensing at the 
outskirts) and of R (e.g. via member galaxy kinematics) 
could constrain the actual effective value of F%, and enable 
extrapolation of the halo radius 7?. 

Figure [5] illustrates the radial structure of our base- 
line minimal-m* solutions with rh = 10 m© yr _1 and 
Tr = 1 keV, but differing in F2. Figure [7] depicts compa- 
rable models with a stronger gas inflow, hotter surface and 
smaller radius. Figure [8] shows solutions with weak inflow, 
and a cool surface at large radius. The gas density profile 
is monotonic in radius, approximating a broken power-law 
with the break appearing at kiloparsec scales, and a slightly 
shallower slope on the outside. In the outer parts, the index 
of w — 1 is co nsistent with the simplest early models of cool- 
ing flows (e.g.lCowie fe Binnevll 19771 ; iFabian fc Nulsenll 19771 : 



iMathews fc Bregmanlll978l ). As expected from analysis, in 
the innermost regions both the gas and halo have singu- 
lar density profiles approximating a Bondi accretion flow, 
Pi oc r~ F '/ 2 for both gas and dark matter. The dark cusp is 
radially smaller than the gas cusp. Note that the dark cusp 
emerges for different reasons than the cusps of hypotheti- 
cal collisionless halos in N-body simulations. The dark cusp 
emerges as a self-consistent, time-independent, hydrostatic 
response to the central mass m, and the gaseous mass inflow. 
In the l atter sense it is related t o an "adiabatic contraction" 
effect, (|Blumenthal et al.lll986T ). It is not a time-dependent 
relic of cosmic structure formation or merger history. 



The dark halo density slope varies more than that of 
gas. A core of approximately constant dark matter density 
surrounds the cusp, spanning from ten to hundreds of kilo- 
parsecs' radius. The halo outskirts are a rapid decline to zero 
density at the surface. As in gasless models ( Appendix (D2|, 
halos with fewer dark degrees of freedom exhibit a larger 
core. The density gradients of the dark matter in the core 
are flattest in the cool, puffy solutions (e.g. Figure [8} than 
in the more compact cases, where the slope is appreciably 
nonzero (e.g. Figure [7}. 

Locally, the dark matter density outweighs the gas in 
some but not all layers of each cluster. Constancy of m at 
the dark matter surface means that pi > P2 in a thin sur- 
face layer. However the halo density dominates gas through- 
out most of the volume, as far inwards as the core radius 
and deeper. This halo-dominated layer is thicker when F2 is 
larger: e.g. reaching in to r ~ 2 kpc in the F2 = 9 model 
shown in Figure [5] For modest F2 (< 8), gas dominates the 
Bondi-like accretion region of the deep interior. For larger 
F2, the steep dark matter density cusp dominates over gas 
in the innermost regions. In the F2 = 9 case shown, we have 
P2 > pi on sub-parsec radii near the origin, beneath a gas- 
rich layer several kiloparsecs thick. The layers dominated 
by dark matter density are smaller for the hot, compact so- 
lutions (top row, Figure [7] than for cool, extended clusters 
(Figure [8]). In principle, this might become observationally 
apparent in cD galaxy kinematics if the gaseous, stellar and 
dark components could be distinguished perfectly. 

The gravitating mass at the outskirts is predominantly 
dark, rri2(r) > mi(r), which follows naturally from the as- 
sumption of cosmic baryon fraction. For models with small 
F2, (e.g. F2 — 3 in Figure|6]), the gas mass dominates within 
r < 100 kpc, and the central mass m, is the dominant com- 
ponent farther in (e.g. r < 10 kpc in this example). Since 
halos with larger F2 are more concentrated (smaller Ri 2 / R) 
the dark-dominated part of the mass profile is thicker for 
greater F2 (middle and right upper panels of Figure [6j . In- 
creasing F2 shrinks the minimal m*, so the central object 
becomes less dynamically significant too. 

Figure [9] shows radial variations of the logarithmic in- 
dex of the total density (pi + p2) for some of the (m*,T)- 
optimal solutions. Cases with lower F2 are flatter out to 
larger radii, as the dark core is larger. The gas inflow tends 
as pi oc r _1 (or steeper) near the centre, and this contribu- 
tion prevents the overall index from reaching zero exactly, 
even deep within the dark core. The steepest index is —4 or 
lower, occurring where the dark fringe drops. Gas dominates 
increasingly at larger radii, bringing the index up to ~ —1.3 
near the halo surface. From panel to panel in Figure|9] curves 
of a given F2 but different [rh, Tr) look alike except for a ra- 
dial dilation. With T implicitly fixed, the halo parameter F2 
controls the proportions of the core relative to the outer sur- 
face R (see radii ratios in Table [1] and Appendices ID 1ID2)) . 
Thus, observing a few signature radii of a real cluster could 
constrain its F2 and R. As a consistency check, satellites or- 
biting beyond R ought to exhibit Keplerian motion. We may 
disfavour some solutions based on their radii: the family of 
(1000,6) models are smaller than 1 Mpc; while the bloated 
family (1,0.4) exceed 10 Mpc (too cosmologically large). 

Many X-ray, kinematic and gravitational lensing obser- 
vations of clusters find total density indices dropping with 
r from tens of kiloparsecs outwards. X-ray analyses typi- 
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cally assume hydrostasis, and probe out to radii of a few 
hundred kiloparsecs (or ~Mpc in rare cases). Strong gravi- 
tational lensing also constrains mass profiles out to 10 2 kpc 
radii, while weak lensing gives statistical evidence at Mpc 
scales. In many instances where measured indices happen 
to range from —1 to —3, this is prese nted as support for 
NFW-like profiles (|Navarro et al. I ll99rj . see Appendix |D7)) . 
Our reference models show similar indices at comparable 
radii, but with crucial differences in the core (natter) and 
fringe (steeper, then the halo truncates). The present paper 
does not attempt to fit specific clusters, but will compare 
model properties to results in observational literature. Such 
comparisons are tentative: the commonly assumed /3-model 
density law (Appendix lD5|l may over-flatten the central gas, 
overestimating the halo density. Likewise, gravitational lens 
models involve subtle degeneracies (e.g. ISaha et al.l 120061 ; 
iLiesenborgs et aT] 120081 . and their references) that might 
confuse cores with cusps. 

Dipolytrope models can naturally describe those clus- 
ters observed to have soft cores. Some X-ray deprojection 
studies have fitted parametric halo models assuming a den- 
sity cusp, and a few appear s ofter than CDM predictions. 
iKatavama fc Havashidal i|2004) find an index 0.47 ± 0.31 in 
the cental 100 kpc of A1835. 



Etto ri et all (|2002al ) indicate 



0.6 in A1795. lEttori et al.l (|2002bT ) prefer a modified Hub- 



ble model (SJD5J) over NF W fits for 10/2 clust ers, giving 
r B of a few hundred kpc. IVoigt fc Fabianl (|2006l) foun d in- 
dices < 1 in 4/12 of their clusters. I Zhang et"aH (|2006l ) fit a 
wide scatter of flattish cusp indices among 13 clusters. Com- 
bining gravitationa l lensing with stel lar kinematics of the cD 
galaxv. lSand et all d200ll2004|2008f) find indices ~ 0.6, and 
< 1 confidently. Rzepecki et alj ( 20071 ) finds an index « 0.7. 
These results at radii < 100 kpc are consistent with the 
shallow regions (r < Rt) of many curves in Figure [9] How- 
ever the most bloated family (upper panel, m = 1 m© yr _1 , 
Tr = 0.4 keV) flattens through Ri at Mpc scales, (implau- 
sibly large) disfavouring cases with F2 ^ 9.5. Among the 
compact family (lower panel) the F2 = 9.5 and 9.9 curves 
are too steep in relevant ranges. The medium cases (middle 
panel) or their homologous relatives (Appendix|B]) are more 
likely representations of normal cored clusters. 

Many gravitational lensing studies treat cluster cores as 
pseudo-isothermal spheres (PIS, Appendix ID4|) . which bear 
comparison to dipolytrope cores. Given any empirical PIS 
core size r s , we can infer R2 > Ri ~ r s . Then one can ex- 
trapolate R from T able [T] ratios, and th is should enclose the 
observable cluster. iDahle et alj (|2003T l found r s ~ 66 kpc 
for an ensemble of clusters. Appraising that this core is too 
small, they rejected fluid-SIDM models of the day. However 
in our calculations this core size predicts a plausible halo 
surface radius of R > 1.9, 3.8 or 17 Mpc (for F 2 = 9.0, 9.5 
and 9. 9 respectively). PIS and NFW fits by lEttori et all 
|2002al ) imply Rt ~ 0.1 Mpc and R 2 ~ 0.49 Mpc in 
A1795. iDiego et al l (|2005T ) fit r s w 15 kpc for A1689. This 
is awkwardly sm a ll, fav ouring higher F2 ~ 10. However 
iBroadhurst et al l (|2005f ) found a core ~ 200 kpc in the 
same cluster, a nd a fringe tru n cating around 2 Mpc (favour- 
ing F2 ~ 7). lHalkola et al.l (|2006l ) found a similar pro- 
file, but with r B ~ 66 kpc, which const rains (F2,R) like 
IDahle et all (|2003l ). iLimousin et all (|2007h fit two halo core 
elements with r s w 99 kp c and 66 kpc . For RX J1347- 
1145, a very massive cluster. lHalkola et alj (|2008h find a core 



r B « 117 kpc while iBradac et all (120081 ) fit 160 kpc. A core 
larger than A1689's befits a heavier system, assum ing uni- 
versal values of F2 and Ri /-R. lRzepecki et all (|2007i ) find flat 
cores in RCS0224-002, with r s « 112 kpc and 12 kpc. We in- 
terpret the larger measurement as the true hal o core, and the 
smalle r feature as baryon-induced contraction. ISaha fc Readl 
(|2008l ) find cuspy profiles for r > 25 kpc in ACO 1703, but 
possibly a density shelf at « 100 kpc. They suggest the shelf 
is meso-structure; we suggest an innate core with a partly 
contracted interior. 

Simple collisionless dark matter models predict p2 oc 
r -3 asymptotically forever, whereas we predict ever steeper 
indices until finite truncation at megaparsec scales. In fu- 
ture, cluster outskirts will become key observational tests 
of halo models. Deeper exposures from newer-generation 
X-ray observatories are needed, and more conclusive grav- 
itational lens models would help. Alrea dy there are lens- 
ing hints of outskirts s t eeper than NFW ([Diego et al.l 120051 ; 
IBroadhurst et alj|2005l;l Umetsu fc Broadhurstll2008l ). An X- 
ray deprojection by iNevalainen et alj ( 19991 ) implies a pe- 
riphery with p2 oc r~ 4 . A much wider cluster X-ray depro- 
jection out to r w 1.7 Mpc (|George et alj|2008l ) shows that a 
hydrostatic NFW-based model cannot fit observations, be- 
cause of an excess of mass or a deficit of gas pressure in the 
fringe. We interpret this as evidence for a finite polytropic 
halo. 



3.3.2 thermal structure 

For dark matter, a\ decreases monotonically with r (e.g. 
dark curves in lower panels of Figures [618) . The temperature 
peaks (like a\ oc r ,_1 ) in the inner density spike, then stays 
flat over several decades in radius, but plummets between 
the core and the dark surface. The upturn from core to spike 
occurs at what is effectively a central gravitational sphere of 
influence, r ~ Gm/a^ Note that the central point-mass m* 
doesn't yet dominate at this radius: the intervening gas and 
halo masses also contribute. Thus the border of the spike is 
essentially a sphere of sei/-influence, where the self-gravity 
terms become important. 

The gas temperature gradient is negative in the cusp. 
Accretion power dominates over cooling, and the tempera- 
ture profile evolves a Bondi-like slope, a\ oc r _1 . Near the 
central object, the ratio of dark to gaseous temperatures 
approaches a constant. Taking the limit r — ► in equations 
(Ml and J39} gives 



2Gm» 



02* 



F2 + 2 

Fi +2 
F2 + 2 



and 



1 + 



Ft 



Fl Ml 



(46) 



(47) 



explicitly involving the central mass and gas Mach number. 

Outside this hot, parsec-scale accretion zone, the gas 
thermal structure depends upon the local balance of radia- 
tive cooling, compression and accretion power. Gas tempera- 
ture gradients may be either positive or negative, depending 
on which terms dominate equation (I20|l . Our optimal solu- 
tions show a local maximum temperature (T max ) at a radius 
typical of the dark core, and a global temperature minimum 
(Tmin) somewhere in the kiloparsec-scale interior. The region 
between these extrema, where dT/dr > 0, is one reasonable 
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Figure 6. Steady inflow structure solutions for models with with the same mass, m = 10 rag yr _1 , Tr = 1 keV but F2 = 3,6,9, in 
columns from left to right respectively. Gas and dark matter properties are marked in grey and black respectively. Top panels show the 
masses interior to radius r, with the central mass m, marked by a dotted line. Middle panels show densities, compared to the cosmic mean 
density (dotted line). The inertial "core" radius Rj and the slope-2 radius R2 of the total density are marked with arrows. The bottom 
panels are temperature profiles. The natural radial unit depends on gravity and bremsstrahlung constants, 1U X = B/G w 0.246 Mpc. 
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Figure 7. Density (top) and temperature (bottom) profiles as in Figure [6] but with rh = 1000 itiq yr 1 and Tr = 6 keV. These are the 
hottest and most compact models. They have the smallest temperature ratio T max /T m i n of the (m„ , T)-optimised models. 
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Table 1. Parameters and global properties of the minimal-m* cluster models. We fix a fiducial total mass, m = A0U m ~ 3.57 X 10 14 mQ, 
and seek cosmic composition, 1/T « 0.163, inside R. From left to right, the columns are: dark degrees of freedom, inflow rate (m© yr~ 1 ); 
surface gas temperature (keV); outer radius (U x ~ 0.246 Mpc units); gas, dark and total concentrations; gravitational radius ratio; four 
of the density slope-radii; the peak of the rotation curve; and the minimal central mass (given in solar units and as a fraction of the 
cluster mass). The signature radii (Ri 1 , Ri 2 , Ri, R w , Ri , R2, R3, R4, Ro) are defined in Appendix IDll 
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definition of the "cooling core" . The peak typically appears 
on the order of a few times 0.1 Mpc, as in X-ray observations 
since the earliest studies of clusters. This characteristic ra- 
dial scale ~ Ux — B /G may be natural to bodies governed 
by gravity and bremsstrahlung radiation (see Appendix IX)) . 
The temperature peak occurs at smaller radii in the com- 
pact solutions than in wide clusters (compare Figures [7] to 
[5| The temperature dip appears at around r ~ 1 — 10 1 kpc, 
with little sensitivity to m. The dip radius is outside both 
the halo density cusp and the break-radius in the gas density 
profile. 

Most of the cooling core overlaps a deep layer where 



dark matter is hotter than gas (a\ > a 2 ). This layer com- 
prises the dark core (except the central spike). Here, any 
disturbance from the cluster equilibrium is likely to cause 
waves of adjustment that propagate faster via the halo than 
gas acoustic modes 

In all our solutions, the gas temperature drops off at 
large radii, as expected in any well bound polytrope. Cooling 
is ineffective at the low densities on the fringe, and at large r 
the velocity terms may also vanish from (|24[) . This leaves the 
gravity term dominant, which guarantees dT/dr < 0. The 
chosen boundary values of Tr correspond roughly to bound 
or accretion- warmed configurations f i]3.2[) . We predict that 
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Figure 8. Cluster profiles as in Figurc[7]but with rh = 1 rriQ yr and Tr = 0.4 keV. These clusters are very radially extended, and 
have a large temperature variation (T max /T m ; n ) in the ICM. 



all isolated clusters have a temperature decli ne at suffi- 
cient l y large radii. Many observations agree (e.g.[Markevitchl 
19981; iFinoguenov et ail l200ll; |Pe Grandi fc Molendil l200a 



Piffaretti et all 120051; IVikhlinin et al.1 120051; IZhang et all 



20061 ; iPratt et alj|2007f). In cases where outer temperature 



profile seems flat (e .g. lAllen et al1l200ll ; iKaastra et"alll2004 
lArnaud et alj|200a ) we predict that a decline will eventually 
appear farther out. 

The finding of T m i n / is a highlight of the model. 
There exists a non-zero floor temperature for every steady 
cluster. Gas does not cool indefinitely, and we have no 
need to invoke distributed mass dropout. It is no surprise 
that observed cooling cluster cores lack massive conden- 
sations of cold gas and extragalactic stars. More signifi- 
cantly, the existence of a temperature floor does not require 
non-gravitational heating. It emerges simply from the co- 
adaptation of gas and halo profiles in their shared gravita- 
tional potential (whether an active galaxy operates or not). 

When constrained to cosmic composition and minimal 
m t , both T max and T m i n increase with increasing F2 (Fig- 
ure [10]). However the ratio T max /T m i n appears less sensitive 
to F2 than to the other global parameters, such as R. In 
the moderate models (with rh = 10 tuq yr -1 and Tr = 1) 
the ratio is 10 < T max /T m i n < 41, and usually « 35. For 
the smaller, strong-inflow models (with rh = 100 ttiq yr -1 , 
T a = 3) we find 8 < T max /T min < 22, For the most com- 
pact and heavily accreting model (with rh = 1000 rriQ yr -1 
and Tr = 6) we find 5.3 < T max /T min < 8. At the oppo- 
site extreme, for the radially largest, coldest series of solu- 
tions with lowest gas influx (1 yr -1 , Tr = 0.4), we find 
a deep temperature contrast: 35 < T max /T m i n < 81. The 
widest clusters provide the greatest cooling length before 
the inflow reaches the accretion-warmed interior. Relatively 
high temperature, compactness and heavy inflow yields the 
smallest temperature variation in the model ICM. 

X-ray cluster observations reveal T max /T m i n ~ 3 to 
4, which is softer than the temperature range of our 
optimal models. However this difference is reconcilable, 



since Tmax/Tmin ratios vary across the solution space. The 
minimal-m» model typically has a temperature ratio near 
the upper extreme. Neighbouring solutions with greater Mr 
or m* can have a warmer thermal dip, and the soft limit of 
7max/Tmi n w 1 is attainable for ?n* several times heavier 
than the minimum. 

One may also ask how far the observed Tmax/Tmin could 
underestimate actual ratios due to finite radial binning and 
imperfect deprojection. The thermal minimum is a thin layer 
at small r; and an observational annulus superimposes hot- 
ter gas from outer shells. If a shell of density pi and temper- 
ature (Ji overlaps the annulus by area 8 A then its emission 
weight is tx p\oi(8A). Figure [TT] shows the contributions of 
shells to two annuli centred on the temperature dip and peak 
of a cluster. For annuli of relative radial thickness ±10%, 
the integrated, weighted temperature ratio (T max )/(T'min) 
shrinks by < 60% (see Figure [121 for examples). Thickening 
the annuli makes little more difference. Projection effects 
cannot wholly hide the strongest thermal contrasts. 

Ratios of Tmax/Tmin ~ 3 may occur naturally if clusters 
have non-minimal m* and the gas physics varies from our 
ideal. Raising the ICM effective heat capacity, Fi > 3 (de- 
scribing micro-scale turbulence or a cosmic-ray component) 
may help. Conduction and non-gravitational heating may 
play a role, though less influentially than popularly thought. 

Bremsstrahlung radiative cooling becomes locally dy- 
namically significant on timescales of 

. _ Pi/(7i - 1) Fid 



2B Pl 



(48) 



In our optimal models, the radial profile of the cooling time 
is approximately a power-law, t c m t £ *(r /r*) Q (see Fig- 
ure [13]). The constant t c * is approximately the same within 
families of clusters with the same (T, m, rh, Tr) but differ- 
ent F-2.. For minimal-m* solutions, the index is fairly con- 
sistent from kiloparsec to megaparsec scales, a m 1.6 — 1.7. 
This is mildly st eeper than X-ray evidence. Fo r example, 
IVoiet fc Fabianl (|2004h and iBauer et all l|2005l ) deproject 
tens of clusters observed out to moderate redshifts, show- 
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Figure 9. Profiles of the radial log-index of the total density 
Pi + P2, for (m*, T)-optimal solutions with F2 values annotated. 
Panels from top to bottom show families (fn/rriQ yr 1 , Tjj/keV) = 
(1,0.4), (10,1) and (100,3) respectively. 



ing 1.3 < a < 1.5 in strata from 10 kpc to 0.5 Mpc. Their 
normalisation also seems comparable to our standard scal- 
ing, since their curves also cross the Hubble time around 
~ 0.1 Mpc. This is a fair agreement, especially considering 
the constancy of the fiducial mass that we imposed, and 
the difference between their assumption of hydrostasis and 
the use of the full Euler equation. Improvement might be 
possible if we were to vary the state of the gas Fi > 3. 

Cooling is cosmologically relevant wherever t c is shorter 
than the Hubble time. This occurs in some sufficiently dense 
inner zone of each cluster. If we define this region's cool- 
ing radius, R c where t c < tn, then we typically find that 
0-2 < R c < 1 (about 50 kpc to 250 kpc, Figure HU). The 
"cool core" can reasonably be defined as the layer r < R c 
rather than in terms of thermal gradients. The cool core is 
usually smaller than the halo core (R c < Ri 2 ). Along the 
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Figure 10. ^-dependency of the peak and dip temper- 
atures of the (m* , T)-optimal models. Black lines show the 
ICM fringe peak temperature. Grey/cyan lines show the inner 
dip temperature. The families of solutions are symbolised by: 
(T R /keV,m/m e yr- 1 ) = (6,1000) (O solid); (3,100) (A dot); 
(1, 10) (□ dash) and (0.4, 1) (Q dot-dash). 



(m*, T)-optimised tracks, we find that R C (T,F2) shrinks 
with increasing F2, despite the increase of R(T, F2) with F2 
and the near constancy of the gas concentration (Ri 1 /R). 

Although formally the cool core is a minor part of the 
cluster volume, cooling controls or affects the exterior gas 
structure indirectly. Contraction and subsidence of cooling 
gas reduces pressure support and draws in the effectively 
adiabatic gas farther out. (Thus the flow does n't need break- 
ing around Rc , as in iBinnev fc Cowiel Il98ll ) . In steady so- 
lutions, the core inflow matches the global rate rh of cos- 
mic accretion from outside the halo. If the central inflow 
caused by cooling doesn't match the global inflow at the 
outer boundary, then a corrective acoustic wave, rarefaction 
wave or shock must propagate outwards into the external 
cosmic medium. Pedantically, it would be misleading to de- 
scribe a cooling flow as driven by external pressure; the core 
slumps because cooling undermines local hydrostasis, and 
the outskirts merely follow in sinkage. 



3.3.3 entropy profiles 

The gas entropy profile is potentially an important diag- 
nostic of the structure and history of a cluster. For a set- 
tled, convectively stable system, dsi/dr > everywhere. 
Hot bubbles float (and cool clumps sink) whilst changing 
volume adiabatically, until settling at a level with compara- 
ble si. Spherical adiabatic accretion is expected to yield a 
power- law slope, si oc r a . Central non-gravitational heating 
may create a radially constant "entropy floor" comprising 
an "isentropic core." For gas that shocks as it accretes into 
a cuspy halo, theory and hydrodynamic simulations predict 
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Figure 11. Any observational projected annulus includes gas of 
temperatures (T) superimposed from a range of 3D shells. We plot 
emission weights of gas shells crossing two mock-observational 
annuli: the grey curve depicts matter in an annulus around the 
thermal dip r £ [0.9, l.l]iJ m i n ; the black curve shows an annulus 
around the thermal peak r £ [0.9, l.l]-R m ax- The model has m* = 
2.85 X 10 7 m Q , F 2 = 9, T R = 1 keV and rh = 10 m Q yr" 1 . 
Projection reduces the apparent peak/dip temperature ratio from 




T max /T min a 5.96 to 3.96. 
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Figure 12. Comparison of actual ICM temperature range 
Tmax/Tmin and the projected ratios (Tmax)/(Tmin) calculated 
with weights as in Figure [TTJ with the same (Fg, m, Tr). Each 
point is a model with non-minimal m* near the "cold" border. 



a ~ 1.1 and no fiat core (|Tozzi fe Norm an 2001; , Kavll2004 
IVoit et alj|2005h . 

Our solutions show different indices in distinct layers. 
In the cooling core (10 kpc < r < R c ) we find a ~ 1.7, 
but a < 1 in the effectively adiabatic outskirts. Our model 
lacks non-gravitational heating, so our solutions never de- 
velop a constant entropy floor. Instead, the entropy profiles 
soften to a ~ 0.2 in the hot accreting interior around/within 
the outermost (quasi-) sonic-point (r ~ 10 kpc; see e.g. the 
bottom panel of Figure I15[l . This slope persists inwards for 
several orders of magnitude in radius. The flatness means 
that cooling is less significant in the warm kpc-scale inte- 
rior than in the cool-core, but there still remains a shallow 
stratification of si. 

The entropy slopes in the ICM are not very different 
from observation, though less agreeably than the indices 
of the cooling time. Flat isentropic cores are not observed 
jEttori et all l2002al : IPonman etakl 120031 ; iPratt fc Arnaudl 



Figure 13. Cooling timescalc profiles of (m* , T)-optimal clus- 
ter solutions, for F2 = 9 halos. Curve patterns correspond to the 
(rh, Tjj) cases in Figurc lTuI The grey horizontal line marks a time 
of 13.7 Gyr, assuming an ion composition factor £ = 1, as in Ap- 
pendix [A] We assume the default cluster mass scale m = 40£/ m . 




4 6 
halo degrees of freedom F z 

Figure 14. Variation with F2 of the radius R £ , within which 
cooling is cosmologically relevant, i.e. the cooling time matches 
the Hubble time. We have assumed the cluster mass m = A0U m - 
The symbols represent the same families of solutions as in Fig- 
ure [10] The hot /compact, large-m clusters have the largest R c . 
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120031 ; IPiffaretti et ahl l2005l l , though there is debate about 
whether the entropy normalisation of smaller systems is af- 
fected by to feeback heating or preheating. The entropy 
ramp appea rs at si ~ 20 — 140 keV cm 2 in observe d 
profiles (e.g iLlovd-Davies et al.1 l200d : iDavid et all l200lf ). 
It begins at ~ 1 keV cm 2 in minimal-m, models, but 
higher for non-minimal cases (e.g. 10 keV cm 2 for that 
in Figure lll[) . Crudely, we expect the stellar matter po- 
tential to have a raising effect similar to large m* (work 
in preparatio n). Most au t hors fi nd slopes a ~ 1 in th e 
cooling core lEttori et al.1 (120023 ) ; IPiffaretti et al.l l|2005l) ; 
IPratt et al.l lj200fj ): IZhang et al.1 l|2006l ) found a « 0.97, 
q ~ 0.95, a ~ 1 .1 and a = 0.99 ± 0.06 respectively. 
iLemze et al.l (|2008i ) combined X-ray and lensing maps to 
model the mass of A1689, assuming spherical hydrosta- 
sis and outskirts declining as pi,p2 oc r~ 3 . In the range 
10 kpc < r < 1 Mpc, they exclude any entropy floor, and 
found indices a = 0.82 ±0.02 when the halo was freely fitted 
(or 1±0.2 when forcing a cuspy halo model). Thus observed 
slopes are slightly flatter than either simulations or our so- 
lutions. This may be an artefact of the NFW-like profiles 
assumed in the analysis of data and construction of simu- 
lations. Alternatively, gas physics with F\ > 3 or thermal 
conduction might improve the match. 

Several X-ray observations have seen a softer entropy 
slope in regions inside r < 20 kpc, and the plots show an in- 
dex 0.1 < a < 0.5, in agreement with the results here. Such 
observations are uncomm on as they requ i re fine resolution of 
thin annuli at small radii. IPonman et al.l (|2003l ) plot entropy 
profiles of 66 varied objects, and the curves seem softer at 
smaller radii. Though unc l ear, t his might be the start of an 
a < 0.5 ramp. IPratt et al.l| 2006h observed mild flattening at 
r < 10 kpc. IDavid et all l)200lf ) deproject the Hydra A clus- 
ter finely in over 30 annuli, showing a < 0.5 cl early across 
the innermost four of them. I Russell et al.l (|2007t ) observed a 
cool-core group with no active AGN, and an entropy slope 
of 0.5 appears within r < 15 kpc. When this kind of entropy 
break appears and attracts comment, AGN heating is con- 
ventionally invoked. The observed shallow region has been 
regarded as merely the outer edge of an (unseen) slope-zero 
isentropic interior. We argue that this slope is actually a 
signature of subsonic inflow inside the (quasi-) sonic point 
(i |3.3.5|l . and not necessarily due to heating. We predict that 
the a « 0.2 zone persists inwards to the nucleus. 

Nevertheless we caution that the inner entropy slopes 
may prove hard to test. Fitting a flat-cored /3-model could 
potentially underestimate the central gas density, since a 
classic cooling flow has pi k r _1 at relevant radii. Tempera- 
ture variations can hide this slope, giving an apparently flat 
X-ray brightness core. Underestimating gas density could 
lead to overestimates of the entropy, exagerrating the flat- 
ness of an observed inner Si(r) profile. 



3.3.4 illusory mass deposition 

Early X-ray imaging studies used approximate formulae to 
estimate m from their deprojected temperature and lumi- 
nosity profiles. In the notation adopted in this work, the 
gradient equation for gas temperature (|20|) was commonly 
reduced to 
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Figure 15. Radial profiles reaching the deep interior, for 
(m», T)-optimal solutions where rh = 10 tuq yr -1 . Tr = 1 keV 
and F2 = 3, 8, 9, 9.9 (annotated). Panels from top to bottom show 
Mach number squared, enclosed non-gaseous mass, and gas en- 
tropy. 



rh(r) 



L«r) 



\a\ + $(r) — $(inner) 



(49) 



where L is the luminosity emitted by a nnuli interior to r. Ki - 
netic terms are dropped. For example, IStewart et al.l (|l984l ) 
plotted m(r) dropping towards the centre. This was taken 
as evidence for mass dropout throughout the cooling core, 
and was justified in terms of runaway local thermal insta- 
bility in a multiphase medium, depen ding on some initial 
spectrum of dumpiness l)Nulsenlll986l ). Radially varying m 
and multiph ase gas became standard ingredients of co oling 
flow fits (e.g. lThomas et al.lll987l:IJohnstone et al.lll992l ) and 
theories (e.g. IWhite fc Sarazi iil ll987l ). However the implied 
deposition products (cold gas and stars) are not observed in 
sufficient amounts, leading to the "cooling flow problem." 

When we apply formula (|49|l to the gas density and tem- 
perature solutions, it reproduces the true m somewhere near 
r < Rc, but at smaller radii the formula underestimates rh 
(Figure [16|). The fictitious radial variation of "rh" resembles 
observationally derived curves. Such profiles are caused by 
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Figure 16. Hydrostatic imaging estimates of rh give a misleading 
appearance of radial variation. The true inflow rate (dotted) is 
constant, rh = 10 t?iq yr —1 . We plot cases with Tr = 1 keV, and 
varied halo types: F2 = 2,6,9.9 (annotated). Curves are plotted 
up to the cooling radius R L . 

omission of velocity terms. It is a signature of proximity to 
the kpc-scale (quasi-) sonic point. Thus we propose to rec- 
oncile centrally depressed "to" profiles with the dearth of 
massive cold condensates on cluster scales. Mass dropout 
from the keV-temperature cluster medium may be unneces- 
sary, and a single-phase ICM may be plausible after all. We 
conclude that that gas inflows reach the central object when 
the system is settled, or perhaps stall at small, intra-galactic 
radii during any episodes when the structure is driven briefly 
off stationarity. We discuss this further in ij4.ll 

We have also tested the effect of the hydrostatic ap- 
proximation on mass profiles estimated from gas density and 
temperature observables. We find that the errors are neglible 
outside > 10 kpc, and the mass underestimate is only a few 
percent at kiloparsec radii. 

3.3.5 Mach number profile & bottlenecks 

The Mach- number profiles of (m t , T)-optimal solutions 
show a variety of features almost anywhere from the cluster 
surface to the smallest calculable radii, (e.g. top panel of 
Figure [To)) . There is at least one local maximum in M 2 , and 
they are more numerous (per decade in r) if F2 > 8. In many 
cases, the flow becomes transonic (M 2 — > 1) sharply around 
one of the maxima. This is corresponds to the "sonic point" 
of a maximal subsonic solution in simpler m odels of adia- 
batic, non-gravitating accretion (|Bondill 19521 ). A slight vari- 
ation of the outer boundary conditions can change a sonic 
point into a supersonic break. Thus these points define bot- 
tlenecks that incur the "too fast" border, (as in Figure [TJ 
thereby constraining the set of steady solutions. The radi- 
ally outermost M 2 maximum is often (but not always) the 
tightest bottleneck. Throughout our survey of rh and Tr pa- 
rameters, the outer bottleneck usually occurs somewhere in 
the radial range 1CT 4 < r < 1(T 2 (~0.25 - 25 kpc). 

Bottlenecks in the gas entropy si can also occur, where 
small variations in outer boundary conditions can trigger 
a cooling catastrophe in a specific radial band. The most 
susceptible radial layers define the "too cold" border. 

The cause of bottlenecks appears in the structure of 
the Mach- number equation, (|37[) . The cooling term Plt c 



vanishes at small radii, so that the geometric and gravi- 
tational terms compete to determine the sign of dM 2 jdl. 
Approaching the origin, there is less interior mass, and in 
some solutions the ratio m//3 ai shrinks enough to guarantee 
that dM 2 /dl < 0. As M 2 increases nearer the origin, the de- 
nominators (1 — M 2 ) shrink, which steepens the gradients of 
gas-related quantities. Often this leads to a runaway descent 
into a supersonic or cold catastrophe. In other conditions 
there is a narrow escape and return to low-.M conditions at 
deeper radii. Even in these "narrow escapes" , the behaviour 
in the deep interior can be dominated by "see-sawing" be- 
tween positive and negative terms in the ODEs. In the set 
of surviving solutions, the multiple "narrow escapes" appear 
as ripples or steps on the density profile. 

For larger F2, the central mass gradients are steeper, 
and prone to yield low m/ f3 ai values, which triggers more 
bottlenecks and see-sawing in layers near the origin. These 
tend to become more numerous and restrictive as F2 in- 
creases. As the middle panel of Figure [T5] shows, when F2 is 
big the halo structure in the nuclear region is a concentric 
set of uniform cores with steep fringes. Undulations appear 
in the gas entropy profile as departures from the typical 
si slope, coinciding with the deep-core mass shells (bottom 
panel of Figure II 5p . 

3.4 Jeans stability 

The outer surface of a finite single-fluid polytrope occurs at 
a radius comparable to the local Jeans radius, 

As equilibrium structures governed by the balance of self- 
gravity and pressure, they are necessarily stable against 
gravitational collapse. Thus these bodies are Jeans masses, 
albeit with non-uniform interiors. 

The gravitational stability of a two-component system 
is not as immediately obvious, with or without complica- 
tions of inflow and cooling. The effective sound speeds and 
densities of the constituents differ at every location, and so 
their Jeans radii differ as well, r n 7^ r J2 . It is possible for 
the dark matter to be locally Jeans-unstable while the gas is 
Jeans-stable, or vice versa. Figure [TH shows the ratio of the 
radial coordinate to the local Jeans scale, throughout a set 
of reference models with cosmic composition. As expected, 
r approaches rj in the fringe, for both gas and dark matter. 
Descending from megaparsec- to parsec-scale radii, each be- 
comes more Jeans-stable. Indefinitely nearer to the origin, 
the gas becomes ever more stable. 

However the central, sub-parsec gravitational stability 
of the halo is more complex and contingent. For low F2, the 
dark halo is Jeans-stable at all radii. In the F2 — 6 case, the 
halo sits at a nearly constant degree of Jeans-stability in the 
sub-parsec interior. For F2 > 6 halo stability lessens nearer 
the origin. For F2 > 8 the halo hovers near marginal Jeans 
stability for several decades in radius near the central mass. 
The upturn towards this condition begins at radii as large as 
10 pc, and is sharper for larger F2. In steady solutions, the 
ratio r/rj within the cusp never exceeds its maximum value 
near the dark surface. Thus the dipolytropic cluster models 
are formally gravitationally stable, and we cannot infer a 
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spontaneous collapse of the dark cusp, without the onset of 
other instabilities. In Sj4T]we will discuss the possibility of 
externally stimulated collapse. 



3.5 Virial scaling 

Our formulation deals with the stationary, relaxed condi- 
tions to which spherical clusters tend ultimately. By con- 
struction, we omit time-dependent, externally driven evolu- 
tion and our results are indifferent to cosmological history. 
Collisionless self-gravitating systems are expected collapse 
to a typical mean density that is some multiple of the cos- 
mic critical density, implying a "virial mass" within a "virial 
radius." The model clusters in this work are not collision- 
less idealisations, but their constituents may have been out 
of acoustic contact before assembly, which is effectively sim- 
ilar. Thus virial scaling may have some back-of-the-envelope 
relevance in observational comparisons. 

Masswise, our solutions rescale freely, and the virial ra- 
dius is not an emergent scale. In order to mimic cosmological 
virial scaling relations, we can rescale the mass, temperature 
and other profiles such that R v matches one of the signature 
scales. In Table we show how three choices of virial radius 
(-Rv = Rs, Rv — Ri, Rv = R) affect the normalisation of 
the (m* , T)-optimal models. We show the mass within the 
virial radius, and the rescaled inflow rate. Our radially ex- 
tended solutions become giant clusters. For them, the choice 
of Rv = R is clearly excessive. The most compact solutions 
rescale to the mass of a group or giant elliptical galaxy. For 
the compact cases with high Fs, the density profile is steep 
almost everywhere, R& is small and thus the choice R v = R3 
implies a puny galaxy mass and negligible inflow. If any sim- 
ple prescription for R v applies to all the model cluster fami- 
lies here, then R v « Ri seems like the most realistic choice. 

Virial scale selection is not necessarily the best choice 
for linking our solutions to a cosmological scenario. We could 
alternatively scale each cluster so that the acoustic crossing 
time is some fraction of the Hubble time. However this range 
of choices is too wide to explore in the present work. For 
the purposes of calculating X-ray spectra (in H3.6l below) we 
shall adopt the simple choice of R v = Ri. 



3.6 X-ray brightness profiles & continuum spectra 

For photons of energy hv, the emission power per unit vol- 
ume due to thermal bremsstrahlung is 



Bpi -hv/tr? 



(51) 



We project the spherical structure of each cluster solution 
by integrating this emissivity along sightlines to produce 
simulated X-ray surface brightness maps. Figure [18] shows 
the brightness profiles of our baseline models in 0.1 — 2.4 keV 
and 2—10 keV bands. The dense central spike appears as a 
bright, steep spot at radial scales r < 10~ 2 U x ~ 3 kpc. This 
is comparable to the size of a cD galaxy, a tiny fraction of the 
cluster volume. The soft X-ray profile declines at all radii, 
but is shallower between the bright spot and R 2 . The hard 
X-ray profile is more clearly structured: steep in the central 
spot; flat in some core ranging from 1 kpc < r < O.I-R2; 
and steep again in the outskirts. It is noteworthy that the 
apparent X-ray core is smaller than the halo core and smaller 



Table 2. Virial masses and inflow rates for some of our (m* , T)- 
optimal solutions, under mass rescaling with three choices of virial 
radius: R v = R3, Rv = Ri and R v = R. Masses are expressed in 
terms of log-^mQ) and inflows in terms of log 10 (rriQ yr -1 ). 
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than the radius where gas temperature peaks. The X-ray 
core contains the radius of the temperature dip. The core 
is smaller when F2 is larger (and seems almost to vanish in 
the F2 = 9.9 solutions). 

Observed X-ray brightness profiles of galaxy clusters 
are empirically fitted with a "/3-model" atmosphere within 
a p resumed isothermal distribution of dark matter or galax- 
ies (|Cavaliere fc Fusco-Femianoll 19761 . see Appendix lD5[) . In 
such fits, the surface brightness varies with projected radius 
b according to 



s(b) = s [i + (b/b c y] 



2,-3/3+1/2 



(52) 



Structure, inflow & central properties of galaxy clusters 19 




10 u 10" l= 10" IU 10~ 3 10" 



Figure 17. The ratio of the radius to the local Jeans radius, for (m* , T)-optimal solutions with m = 10 rriQ yr -1 , Tr = 1 keV and F2 
dark degrees of freedom. Halo and gas profiles are marked in black and grey respectively. A dotted line marks the instability threshold, 
r = r j . In the F2 > 6 cases the arrows mark where the marginally stable inner halo encloses a dark mass of 10 7 , 10 s , 10 9 »tiq . A sufficiently 
large perturbation might make this mass collapsible. 
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Figure 19. Luminosity / temperature relation of model clusters. 
The mean temperature derives from the total mass and internal 
energy of gas (fcT x = |C/i/mi) in the ICM (r > 0.01U X ). Circle 
data are the models from Table [T] Squares, diamonds and stars 
are rescaled to "virial" density according to the rules in Table [2] 



Figure 18. Brightness profiles of X-ray continuum projected 
from (m*, T)-optimal solutions with fn = 10 m© yr — 1 and Tr = 
1 keV, but various F2. Black lines show 0.1 — 2.4 keV emission; 
grey shows the 2 — 10 keV band. We have rescaled the mass and 
temperature to a virial selection of R v = Rj. The X-ray core is 
smaller than the halo core, and both shrink as F2 rises. 



The index j3 fits the outer slope, while the parameter b c mea- 
sures the core radius. The observed range is 0.4 < f3 < 1.4 
and t y pically 0.1 Mpc < b c < 0.5 Mpc (e.g. I Jones &: Formanl 
19841: iNeumann fc Arnaudl Il999l ; lEttori fc Fabianl Il999l ; 
Mohr et al.lll999T l. Variant models exist to fit more centrally 
peaked clusters , by combining c usps and/or (5 compone nts of 
different sizes (|Xu et all 1 199S1 ; iPratt fc Arnaudl |2002| ). The 
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/3-model was originally derived for isothermal atmospheres, 
but it is still applied to clusters that are now known to have 
radial temperature variations. 

Excepting the central spike, the /3-model resembles 
the synthetic X-ray profiles presented here. In the rh — 
10 m© yr _1 series, we find core sizes 6 C < 0.3 Mpc, shrinking 
slightly with rising F%. However the (m», Y)-optimal model 
for F2 = 9.9 forms no X-ray core at all, which is empiri- 
cally unfavourable. The fringe brightness declines with sim- 
ilar slopes in all bands: /3 ~ 1.0 ± 0.2 around R2 < r < Ri. 
These results sit within the observed range. 

We note that /3-like profiles are not a distinguish- 
ing feature of our formulation, nor cored halos as a class. 
Even in the entirely cuspy profiles predicted for collision- 
less dark matter, the gas settles as a cored atmosphere, 
with s hallow central grad i ents, fitting so m e kind of (3- 
model (|Navarro et al.lll996l ; fake et all 119981 '). iMakino et all 
(| 19981 ) have shown that /3-like profiles are broadly natural 
to isothermal hydrostatic clusters, though cuspy halos yield 
smaller 6 C . 

Like the projected light distributions, the three- 
dimensional luminosity profile L = L(r) has a flat core, a 
declining fringe and a bright nuclear spot. The central X-ray 
spike fits inside radii r < 10~ 2 , so we choose this as a cutoff 
defining the luminosity of the intracluster medium. At stan- 
dard scaling, m(R) = 40C/ m , the ICM luminosity increases 
with rh and with F2. We find Li cm > 2 x 10 41 erg s _1 in the 
wide/cool family of clusters (rh = 1 m© yr _1 , Tr = 0.4 keV). 
For the medium-sized solutions (rh = 10 m© yr _1 , Tr = 
1 keV), we find Licm > 5 x 10 42 erg s _1 . In the more 
compact solutions we have Lic m > 1.4 x 10 44 erg s" 1 (for 
rh = 100 m© yr _1 ,TR = 3 keV); and L icm > 2x 10 45 erg s" 1 
for the most compact set (rh — 1000 m© yr~ , Tr = 6 keV). 

The nuclear spot is more luminous, but this depends 
on where we truncate the model. If we had undertaken a 
relativistic formulation f i]4.2.6[) . then structures at the clas- 
sical inner radius r» = would shift out to an event horizon 
at r. = 2Gm,/c 2 . If we truncate L(r) at r ~ 5r. (suiting 
accretion onto a black hole) then we find typical nuclear lu- 
minosities around L„ ~ 5 x 10 44 (m/l m© yr _1 ) erg s , 
which is plausible as the total accretion power of an ac- 
tive galactic nucleus. Realistically, we would expect a large 
share of this power to emerge in forms other than observ- 
able bremsstrahlung radiation. Extra physics and sub-parsec 
AGN anatomy may reduce the emergent luminosity to a 
small fraction of the accretion power (see H4.ll £ f4~2.7|) . In 
particular, Compton scattering and dense obscuring inter- 
stellar media must alter or reprocess the escaping radiation. 
Thus, the central density and accretion power in the present, 
undetailed model may be consistent with radiatively ineffi- 
cient black hole feeding or with AGN, which is not such a 
dire outcome as early reviews assumed. (Indeed the bright 
spot is hot and therefore qualitatively different from the cold 
catastrophe accumulations that early works predicted.) 

All of the ICM and nuclear luminosities given above 
have assumed the standard normalisation, with total clus- 
ter mass m(R) = 40C/ m « 3.57 x 10 14 m©. If we rescale the 
masses by a factor X then the luminosity changes by a fac- 
tor X 5//2 (see Appendix [B]). For instance, if we lighten the 
compact, megaparsec-scale, rh = 1000 m© yr _1 solutions 
by a factor X = 0.01 to represent an isolated giant ellipti- 
cal protogalaxy, then the inflow shrinks to rh = 1 m© yr _1 , 



erg s and ICM lu- 
. The central mass 



with nuclear luminosity L„ ~ 5 x 10 42 
minosity L icm » 2 x 10 40 - 2 x 10 41 erg s 
limit rescales as well: the central black hole must weigh at 
least 1.2 x 10 5 m© < m» < 3 x 10 9 m© (depending on L2). 

Luminosity and mean ICM temperature correlate in our 
basic models, Li cm oc T" (see Figure I19[l . For raw models 
set via inflow condition (145|) . the slope is a_ = 3.46, which is 
steeper than observed (~ 2.6 to 3.0, e.i 
Arnaud fc Evrar d) 1 19991 : iNovicki et al.l _ 
20021 : iLumb et al.ll2004h . Clusters rescaled to virial density 



( MM show a = 2.22,2.02 and 2.05 (for R v = R 3 ,Ri,R 
respectively) resembling t he a = 2 pre diction from grav- 
itational collapse theory ijKaiser! 1 19861 ). Nature's way of 
breaking the mass homology seems intermed iate: the cause 
may i nvolve the boundary condition, heating (|Ponman et al.l 
1 19991 : iLoewensteir] l200oh or other aspects of gas physics 

(EM- 

We have calculated synthetic X-ray spectra at several 
projected radii in the bloated, medium and small cluster 
solutions (see Figure [20}. As with the projected bright- 
ness profiles, each of these plots assumes the virial mass 
scaling R v = Ri. This makes the large clusters hotter 
and more massive than under the original normalisation 
(m(R) = 40[/m) becoming a massive supercluster with tem- 
peratures of tens of keV around the virial radius, which is 
several megaparsecs. This rescaling does not resemble any 
known realistic object; the choice of R v — Ri appears un- 
suited to this family of solutions. Clusters of medium ra- 
dius (middle row) show temperatures of several keV near 
the virial radius, which is near or outside the temperature 
peak. At the projected radius of the temperature dip, the 
continuum curves in a clearly multi-temperatured way. 

The results are similar for the smaller clusters (bot- 
tom row) but the peak temperature is w 1 keV and the 
continua are more obviously curved in the 0.1-10 keV band 
that we display. However line cooling is significant compared 
to bremsstrahlung in low-mass systems such as these, which 
may alter the profiles appreciably f H4.2.1jl . Further calcu- 
lations are needed, with a more detailed cooling function, 
to model galaxy/group spectra including line and edge fea- 
tures. However line cooling would break the homology of the 
present solutions (Appendix [B]), requiring a more expensive 
exploration of parameter-space. We defer this topic for fu- 
ture investigation. 



3.7 projected mass and gravitational lensing 
profiles 

To help predict and assess the gravitational lensing signa- 
tures of our cluster models, we calculate mass maps by inte- 
grating pi + p2 along lines of sight at various projected radii 
b. Figure [21] presents the projected profiles of total mass 
column density, E = £(&), for several represenative (m», T)- 
optimised models. The radial gradients are steep on scales 
from several hundreds of kpc to several Mpc (beyond the 
effective core radius Ri). The distribution flattens on scales 
between tens to hundreds of kiloparsecs. (This is the mass 
core.) In the deepest inner regions, within tens of parsecs 
of the origin, the central spike dominates and the flat core 
steepens again. The border between core and spike typically 
occurs around 10 _4 l7a; — lO -3 ^, regardless of whether the 
cluster is cool/bloated, moderate or hot /compact. The ra- 
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Figure 20. Projected X-ray spectra (continuum only) from the models with fiducial parameters (m/m© yr — 1 , Tji/keV) = (1,0.4), 
(10, 1) and (100, 3) (from top to bottom rows respectively). Here, however, the masses and temperatures are rescaled to virial conditions, 
choosing R v = Rj, as in Table[2] The halo freedom is F2 = 3, 6, 9 in the left, middle and right columns. In each panel, the blue (dotted) 
spectrum reaches the temperature dip radius; red (dashed) reaches the temperature peak; and the black spectra are projected at fractions 
of the virial radius, b = {^ij^i^i^i^il} Rv The bold curve is at b = R v . Rescaled values of the peak and dip temperatures (keV) 
are noted in the top-right corner of each panel. 



dius of the projected mass spike is comparable to, or smaller 
than, the radius where the X-ray brightness spike begins. 
The spike in E exists in theory but may be unobservable in 
practice: dominated by the stellar mass of a cD galaxy. 



As was apparent in three-dimensional density profiles, 
the projected E core is smaller and denser when F2 is 
greater. The outer radius R varies more slowly with F2, so 
the mean slope of the outskirts is shallower if F2 is high. 
The spike profile depends on F2 also: for F2 < 6 the gas 
density dominates at ~ 10 pc scales, and the profile has a 
logarithmic slope of ~ — 1. For F2 ^ 6 the spike is more 
often dominated by dark mass, and E shows a slope steeper 
than —1. These predominantly dark spikes are more preva- 
lent and more radially extended in the compact/hot cluster 
models. In the coolest, most diffuse models (top left of Fig- 
ure l2"Tj) all the spikes are gas-dominated down to 1Q~ 6 U X . 



We predict that gravitational lensing measurements 
that probe images on Mpc scales will see steep mass slopes of 
the halo fringe. Medium-separation lensing systems (probing 
just within the core radius) may reveal shallower gradients 
of the dark core. However E indices are unlikely to reach 
zero in the core, due to the pi oc r^ 1 gas contribution, and 
the projection of outer layers. If the innermost regions can 
ever be discerned through the light of the cD galaxy then 
they may appear spiky. 



4 DISCUSSION 



4.1 implications in galaxy evolution 



As we report in £13.3.51 the bottlenecks that select the valid 
stationary solutions occur on two different radial scales. The 
outer critical region occurs at kiloparsec radii. It may not 
be a coincidence that this is the characteristic extent of the 
stellar matter in an elliptical galaxy. If a primitive, initially 
starless proto-cluster or proto-group were perturbed from a 
stationary state, then the cooling cata strophe will tend to 
eme rge within < 3 kpc f rom the centre. iNulsen et al.l (|l984l ) 
and iFabian et all l| 19841 ) anticipated this size from back-of- 
the-envelope arguments; our radially complete calculations 
substantiate it. During a strong enough disturbance, ex- 
cess cold gas may drop out as star formation, until a new 
stationary configuration settles. In short, bumping or shak- 
ing a proto-galactic globule may spawn a spheroidal galaxy 
monolithically. However the bottleneck radii limit the ex- 
tent of stellar condensation, influencing the high-mass cutoff 
of g alaxies (perhaps a l ongside AGN or conduc t ion effects , 
e.g. ISilk fc Reel Il998l: iBegelman fc Nathl 120051 ; Isilkl |2005l : 
iBest et al.ll2006l : IFabian et al.ll2002t ). 

The inner critical regions occur on sub-parsec scales, 
comparable to the sizes of galactic nuclei. Here again, tran- 
sient disturbances towards overcooling might form a local 
concentrations of stars. However in high-F2 cases, the inner 
halo sits on the brink of Jeans instability ( £13.4|) . and the 
upturn towards this threshold begins as far out as ~ 10 pc. 
In these conditions, a non-linear perturbation may trigger 
collective gravitational collapse by the dark matter (until a 
new stationary solution emerges). Supermassive black holes 
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Figure 21. Projected profiles of the total projected density, £(&), of (m* , Y)=optimal cluster solutions with different F2 values 
(annotated circles). We scale each model to the default total mass, m(R) = 40t/ m . 



may be the natural product. The spikes of high-.Fb halos 
may collapse themselves into holes whenever anything rat- 
tles the central galaxy appreciably. The collapsible spike 
mass is comparable to real SMBH masses. 

Earlier works have proposed mechanisms for black 
holes to feed and grow from self-interacting dark mat- 
ter, evading Eddingt on limi t s and emitting l it tle di rectly 
luminous evidence (ISoltanl Il982f ). lOstrikerl (|2000t ) and 
iHennawi fc Ostrikerl (|2002l ) considered dark Bondi accre- 
tion onto stellar black hole seeds, followed by slower, dif- 
fusive loss-cone refilling. Their feeding scenario requires 
SIDM interactio ns to be weak enough to provid e an initial 
NFW-like cusp. iMunvaneza fc Biermannl (|2005l ) show that 
a SMBH can grow at an asymptotic (exceeding exponen- 
tial) rate if the dark hal o has a degenerate fermion core. 
iBalberg fc Shapiro! (|2002l ) found that gravothermal collapse 
in conductive, weakly-SIDM can spawn 10 6 m© black holes 
directly (and these require gas accretion to reach 10 9 tjiq). 
Like the latter mechanism, our scenario requires no seed, and 
its preconditions are self-organised by the gas/halo interplay. 
Like the fermion-feeding model, we are free to consider in- 
teractivity at a strength that precludes NFW cuspiness in 
any epoch. The mechanism is a differentiated Jeans collapse: 
when stimulated, the hole feeds from a single dark gulp, and 
this intake is only limited by the mass of the dark spike. 

Stationarity of the whole cluster demands the de- 
velopment of a central mass. Note however that the 



m* limits express only the minimal object. This lower 
limit does not predict the correlation between supermas- 
sive black holes and stellar bulge properties. Such rela- 
tions probably involve extra regulatory processes operat- 
ing on galaxy scales. The depth of the galaxy potential 
and the kick velocities of black hole mergers may influ- 
ence t he upper limits (e.g. lHaimanll2004l;lMadau fc Quataerj 
2004 iGonzalez et alj|2007l; ISchnittman fc Buonannol 120071 



Campanelli et al.ll2007a lbl; lBogdanovic et al.ll2007 ; Volonteril 
20071 ). In our scenario, massive black holes can form rapidly 
and darkly. If a merging galaxy is apt to expel its SMBH 
then we would expect a replacement to condense when the 
scene settles to a resurgent gas inflow and dark spike. 

During the tranquil periods between mergers, there are 
two plausible fates for the gas inflow after it penetrates be- 
low kiloparsec radii. It may suffer cold catastrophe (if the 
configuration is near the critical border) depositing cold gas 
and new stars centrally. This essentially shrinks the "cool- 
ing flow problem" down to parsec scales, which is arguably 
an improvement over the old 0.1 Mpc-sized problem. More 
likely, the inflow feeds the central black hole directly. 

The existence o f powerful quasars in the era z > 6 
(|Willott et al.ll2003l ; IWalter et al.l [20031 ; iBarth et alj|2003l ) 

requires that black holes grew rapidly. Our model not only 
allows this, but requires it to happen before a galaxy or 
cluster achieves stationarity. Whether large rh inflows are 
sustainable into the modern era is another matter. The ac- 
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cumulated stellar mass distribution may eventually alter the 
central potential enough to change the allowed domain of 
steady feeding solutions fi )4.2.4[) . If high rh inflows can per- 
sist, the implied AGN activity needs explaining. Much of 
the time, an AGN may accrete in a radiatively inefficient 
mode: many complicating factors could reduce the settled 
AGN luminosities to a small fraction of the actual accre- 
tion power. Opacity (i ]4.2,6|l and thermal conduction ( £|4.2,2[) 
may smudge or dim the nuclear bright-spot. An opaque, 
super-Eddington inflow entrains and s wallows much of its 
own luminosity (|Begelman|[l978l . Il979h . The non-luminous 
power of jets and bulk outflows from sub-parsec regions may - 
hide much of the accretion power (e.g. lDi Matteo et al.1120031 ; 
lAIlen et alj|2006l ). 

If cD galaxies form from massive cooling inflows, then 
they may differ qualitatively from elliptical galaxies in the 
field. We expect an association between cooling cores and 
the presence of a cD galaxy. A cluster merger might displace 
an old cD galaxy, and temporarily disrupt the cooling core. 
A new cD galaxy would eventually sprout in the middle of 
the resuming cooling flow. 



4.2 extra baryonic and radiative physics 

Here we briefly discuss several baryonic phenomena that 
could potentially modify our quantitative results. 



4-2.1 line cooling 

The presence of line cooling would break the 
mass/temperature homology of the solutions assuming 
bremsstrahlung cooling alone. Additional parameters are 
therefore needed in the hydrodynamic formulation, thus 
increasing the dimensionality of the problem and the 
complexity and computational cost in searching for the 
structure solutions. 

Line emission must affect cooler systems (galaxies and 
groups) more than rich clusters. The extra cooling may 
tighten the constraints of the "cold" border, perhaps raising 
the lower limits on m»/m for ligher objects. This may also 
alter the compactness of clusters of given mass, comosition 
and m: i.e. the T-tracks may shift in the (F2,R) plane. 

Hydrogen and helium cooling must have been signifi- 
cant for primordial "mini-halos" - hypothetical, metal-free 
bodies of gas and dark matter that were lighter than modern 
groups or galaxies. If their radii are smaller than the most 
compact models in this paper, then mini-halos develop pro- 
portionally larger central masses. The contribution of line 
cooling above bremsstrahlung may raise the minimal level of 
m*/m further. This supports the notion that massive black 
holes condensed directly from primordial envelopes, and that 
condensation may have been as fast as free-fall or the onset 
of cooling catastrophe. 



4-2.2 thermal conduction 

The role of thermal conduction in th e intracluster medium 
has b ee n contentious for d ec ades ( e.g. [Mathews fc Bregmanl 



19891; iTribbld Il989l: ISuginohara fc Ostrikeil Il99 



Naravan fc Medvedevl l200ll; lLoebl |2002|; IVoigt et all 
20021; IZakamska fc Naravanl 120031; iKim fc Naravanl 12003^ 



Voigt fc Fabianl |2004| ; IConrov fc Qstrikerl 120081 ). Locally 



19781; [Biimev fc Cowiej Il98ll; iNulsen etafl 



Tucker fc Rosnerj Il983l; IStewart etail Il984l; [Friacal Il98rj; 
Bertschinger fc Meiksinl 19861 ; Boehringer fc Fabianl 



1982; 



tangled magnetic fields may hinder conduction by inhibit- 
ing diffusion transverse to field lines. However, fields that 
have been aligned radially by a bulk inflow might rather 
promote radial h eat c onduction |Bregman fc Davidl 1 19881 ; 
ISoker fc Sarazinl Il990t ). Unsuppressed conduction could 
warm the coldest layers and lessen radial temperature 
variations: raising the dip temperature and reducing the 

Tmax/Tmin ratio. 

In the conventional models, cluster outskirts have been 
seen as a plausible reservoir for conductively heating the cool 
core. This was a natural proposition when solutions were 
sought in the "too cold" domain, and the whole core was 
thought to suffer a multi-phase, distributed cooling catastro- 
phe. We emphasise that heat conduction outwards from the 
hot nucleus may also be relevant. These regions are ideal for 
conduction, since the densities, temperatures and tempera- 
ture gradients are all high. Conductive cooling of the nucleus 
must flatten the central temperature and density profiles, 
dimming the central luminous spot. Some of the central ac- 
cretion power conducts outwards, ultimately to emerge as 
ICM luminosity elsewere. Conductive heating may loosen 
the constraint of the "too cold" border, enabling smaller m* 
values. We expect this quantitatively, but it requires proof 
from extended analyses. 

4-2.3 thermal instability 

The relevance of thermal instability to the phase structure 
of cool cores h as also been debated since the early theories of 
coolin g flows (|Fabian fc Nulsen1ll977l ; iMathews fc Bregmanl 
Il978l ). In uniform, isobaric conditions, overdense clumps 
may overcool radiatively in a runaway mann er, and con - 
dense as a cold phase in thin, hot surroundings (|Fieldlll965l ). 
However, dense blobs may rain towards the gravitational 
centre, resulting in ablation and warming that counteracts 
instability. Blobs that fall deep and fast enough may reheat 
via shocks (|Cowie et al.lll98d ; lNulsen et al.lll984T l. Buoyancy 
may inhibit instability from the start, by shifting cool or 
warm blobs to strata of matching entropy. Thermal conduc- 
tion hinders thermal instability by warming nascent blobs. 
Local magnetic fields that enwrap blobs might help isolate 
them, aiding thermal instability and mass dropout. Fields 
that thread blobs and their su r roundings might bin d the 
phases to co-move (|Nulsenll 19861 ; iThomas et atil 19871 ) . 

Studies of thermal instability in interstellar shocks and 
stellar accretion suggest that local, homogeneous, isobaric 
analyses are incomplete or too simple for many applica- 
tions. The shape of the cooling function affects instabil- 
ity; if line-cooling complicates the cooling law then gas 
is stable at some temperature. Macroscopic flow geome- 
try can be influential: in general the Eulerian operator 
(dt + v • V) coup les temporal evolution with motion and 
gradients. (See e.g. Chevalier fc Imamurall982l ; lBertschingerl 
Il986l ; ISaxton fc Wulll999l . and references therein.) A bet- 
ter treatment of thermal instability in clusters requires 
analysis of global st r ucture as well as loc al physics (e.g. 
iMalagoli et aLI 1 19871 ; iBalbus fc Sokerl Il989l ). By retaining 
temporal terms in the governing equations, our model is 
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ready for stability analyses that relate regional thermal 
(in)stability to deformations in the halo and gas profiles (in 
preparation). 

However the strongest original argument for thermal in- 
stability fades in our present results. A full treatment of the 
conservation laws and flow velocity seems to dispel an old 
illusion of radially varying m. Constant-rh models resemble 
observations as well, and so we need not invoke multiphase 
effects, nor mass dropout across the core. If mass drops out 
anywhere (for reasons beyond our model) then it happens 
deep inside the cD galaxy, and we have at least reduced the 
"coo ling flow problem" to a sub-kiloparsec AGN problem 
(e.g. iTabor fc Birmevlll993l ). Without dropout, the cooling 
flow leads into hot spherical accretion feeding the nucleus. 

Having a warm T m i n , our solutions are consis- 
tent with the paucity of cold gas and star forma- 
tion measured at l ong wavelengths (see review by 
iDonahue fc Voitl |2004 ). and also with the X-ray spec- 
troscopic evidence that disfavours widespread mass- 
dropout to very low temperatures j Kaastra ct al.l l200ll ; 
iTamura et all l200ll ; iPeterson et~aT1 l200ll . 120031 ). Observa- 
tions both inside and outside the cooling radius are con- 
sistent with si ngle-phase flows but spatially varying tem- 
perature (e.g. | Bohringer et al.l l200ll; iMolendi fc Pizzolato 



200 ll: 



David et al. 2001; Matsushita et al. 2002: Ettorietal 



2002al ). Thus we justify omitting micro-scale thermal insta- 



bility and retain a smooth, single-phased flow model. 



4-2-4 stellar material 

Our present analyses omit the effects of collisionless, stellar 
mass. This is justifiable on cluster scales, where stars are 
effectively passive tracers, and gas comprises most of the 
baryons. Our formulation also applies directly to primordial 
structures, if protogalaxies detached from the background 
before stars became abundant. 

Stellar mass dominates within the effective radii of mod- 
ern elliptical galaxies, such as the brightest central galax- 
ies grown in cluster cores. A stellar mass distribution (see 
Appendix ID10|) may alter our solutions within the central 
kiloparsecs, by deepening the potential there. The extra ac- 
cretion warming may loosen the "cold" constraints on m*. 
It may also soften the ICM temperature dip. We defer the 
evaluation of gas/dark dipolytropes in elliptical galaxy po- 
tentials for future study. 

If the central galaxy is quiesent then direct stellar in- 
teractions with the inflow are negligible. It was recognised 
early that stellar mass loss in modern elliptic al galaxies is 
weak compared to a cooling cluster's inflow fe.g lNulsen et al] 
1 1984 ) ■ and likely to be smothered in terms of mass, momen- 
tum and energy. In our solutions, the central gas is near local 
virial temperature, already similar to the stellar velocity dis- 
persion. Drag from stellar motions may stir gas locally, but 
the cross-section and covering factor of stars is too small for 
much global effect. Thus the main effect of stars is gravita- 
tional and limited to the centre (unless a starburst erupts). 
Starbursts may affect the early life of a cD galaxy, heating 
nearby gas as an AGN might. 



4-2.5 heating by active galaxies 

We omit AGN interactions from our model, in order to focus 
on the undriven, natural tendency of clusters. Feedback ef- 
fects require extra semi-empirical terms, with a diabolically 
tempting number of adjustable parameters. If jet activity 
or other phenomena drive the cluster medium to convection 
or turbulence, with eddy kinetic energies comparable to the 
gaseous internal energy, then the gas gains extra effective 
degrees of freedom, Fi > 3. This may steepen gas density 
profiles, and shift the "cold" and "fast" constraints on m*. 

Some authors propose that AGN output (such 
as the mechanical power of radio lobe s) can suppress 
or even regulate cooling flo ws (e .g. ITabor fc Binnev 



2001 


; Briiggen & Kaiser 2001 


; Quilis et al. 


2001; 


Ruszkowski & Begelmanl 12002); 


Kaiser & Binnevl 


2003; 



Ruszkowski et al.l 120041 ). There seems enough raw available 
power for low and m edium-mass clusters (jBirzan et al.l 
|2004 I Allen et al.ll2006h . but the sufficiency of AGN warm- 
ing in massive cooling flows is debatable. For AGN to be 
effective heaters, radio lobes must mix with the ICM, or 
the gas must be viscous enough to dissipate disturbances. 
The distribution of AGN power is questionable in some 
systems: strong jets may cut themselves channels out of 
the core, dumping th e ir power ineffectually in the outskirts 
(e.g. iBest et~al"1l2006l ; IVernaleo fc Reynolds! 120061 ). In some 
clusters, ICM metallicity gradients imply that mixing has 
not been thorough. Some cool cored clusters lack nuclear 
activity altogether. 

Because of their low entropy, cooling flows may simply 
sink and slip around the sides least disturbed by AGN chan- 
nels. If so then active cooling clusters may still resemble our 
solutions overall. However if the AGN blocks the inflow near 
the centre then the model needs modification there: an ad- 
ditional distributed heating function in C, and a local drop 
in m. Starvation on sub-kiloparsec scales may dim the lu- 
minous nucleus. The outskirts (beyond the reach of radio 
bubbling) should match the standard profiles obtained here. 

However our adaptive model improves or solves key as- 
pects of the "cooling flow problem" without resort to ex- 
trinsic heating. Considering a responsive halo, and all rele- 
vant momentum and self-gravity terms, we find that steady 
clusters always develop a floor temperature due to purely 
gravitational self-warming. This level is high enough to ex- 
plain the rarity of cold condensates. Consequently, we do not 
need non-gravitational heating to dominate on megaparsec 
scales. We dispel the need for finely tuned, two-way feedback 
between the cooling flow and the heating processes. Stellar 
and AGN heating are demoted to an incidental, intermittent 
role. Heating need not be persistent nor stable. 



4-2-6 relativity, opacity and cosmic rays 

In a relativistic formulation, the inner boundary would shift 
from r, = to the Schwarzschild radius of the central mass. 
The bremsstrahlung cooling function acquires relativistic 
corrections. Gravitational redshift dims the inner boundary. 
Inflow velocities will be naturally subluminal. As the gas 
becomes relativistic in the nuclear regions, it gains more ef- 
fective degrees of freedom, Fi — * 6. This implies a higher 
combined heat capacity than a normal F\ — 3 medium. The 



Structure, inflow & central properties of galaxy clusters 25 



region where the dark spike verges on Jeans-instability may 
enlarge. Otherwise the issues of pressure support, the sub- 
sonic constraint and the avoidance of cooling catastrophes 
remain qualitatively unchanged in a relativistic framework. 

If gas in the central regions becomes Compton-thick 
or otherwise dominated by radiation pressure then this will 
also result in Fi — > 6. This will steepen the gas density spike 
and may promote collapse of the dark spike. The opaque in- 
ner inflow must be radiatively inefficient, enlarging the do- 
main of effectively adiabatic Bondi-like behaviour. We can- 
not presently say in which direction the m* limits change. 
Sub-parsec AGN anatomy may complicate the issue. 

If cosmic rays contribute significant pressure to the in- 
tercluster medium, then they would deserve incorporation 
as a third fluid in our model. Extra source terms enter the 
momentum and energy equations to express cosmic ray dif- 
fusion and heat exchanges with the coterminous thermal gas. 
As a relativistic plasma, cosmic rays have F3 m 6 degrees of 
freedom. The consquences for the innermost structures and 
m* limits may follow cases of a halo with F2 ^ 6. 



4-2.7 angular momentum 

If some of the cluster gas possesses significant angular mo- 
mentum then the inflow could deviate from radial streaming 
at small radii. Some fraction would accrete ont o a small disc, 
of a s ize determined by rotational support (e.g. lNulsen et all 
|jL984j). This may alter the minimal-m, limits slightly, and 
the steady solutions might require some minimal disc sur- 
face density profile. If the disc is viscous then it is only a 
temporary residence for inflowing gas; it ultimately feeds the 
central mass. A disc that accretes enough mass could self- 
gravitate and develop density waves, or fragment. This is a 
recipe for forming spiral galaxies. 

However discs (and filaments) are only viable as long 
as they avoid contact with similar bodies. Collisions and 
stirring by asymmetric substructures can restore the spher- 
ical approximation in the long run, and on the scales rel- 
evant to our model. Turbulence or convection in the inner 
regions could easily destroy or preempt a disc, easing angular 
momentrum efflux and quasi-spherical mass influx. Empiri- 
cially, our spherical approximation is valid as far as relaxed 
clusters and central galaxies are actually roundish. 

A spherical model cannot address the topic of halo 
substructure directly. Collisionless dark matter cosmogonies 
overpredict the abundance of satellites a t galaxy and group 
scales by at least an orde r of magnitude (|Klypin et al.lll999l : 
iMoore et all Il999l . l2000l : iD'Onghia fc Lake! hoO^ rSeveral 
authors have debated whether a (weakly) SIDM halo can be 
lumpy and aspherical e nough as a gravitational lens and host 
of sat el lite galaxies (e.g.pVIeneghetti et a l. 2001; Moore et al] 



20001; iGnedin fc Ostriker! 120011 ; iNataraian et al] |20o3 



Furlanetto fc Loebll2002r i7 Dwarf halos might ablate at dark 



bow-shocks, or (we suggest) might persist as self-bound, 
dark eddies that roll as they orbit in a turbulent background. 
Within a cluster, each galaxy perturbs the cluster profile lo- 
cally, and each subhalo centre probably develops a miniature 
analogue of the spherical inflow solutions. 



4.3 nature of dark matter 

In detail, our solutions depend on the assumption of a 
polytropic equation of state for the dark matter. This 
condition is the emergent equilibrium if Tsallis' thermo- 
statistics govern the halo, even if dark matter lacks non- 
gravitational interactions entirely. If dark matter is a Bose- 
Einsten conde nsate, the n it is effectively polytropic in the 
classical limit (|Sinlll994l ; fGoodmaij l200ch ■ If it is a degen- 
erate fermion gas then it has a more complex equation of 
state, obtainable by integrating a loca l distribution function 
|Munvaneza fc Biermann 20051. 20061). or els e a polynomial 
approximation ( Nakaiima fc Morik awa 20071 ). If dark mat- 
ter is collisionless but ruled by long-range dark forces then a 
more complicated treatment becomes necessary, analogous 
to collisionless plasma physics. If dark matter feels strong 
enough local self-interactions, then it is analogous to an adi- 
abatic ideal gas, and a polytropic equation is expected. 

The possibility of loc al dark s elf-in te ractivity 

llSpergel fc Steinhardtl l200d : iGoodmanl l2000l : IPeeblesI 
|2000| ) is theoretically and observationally attractive, but 
not yet exhaustively tested. This scenario explains the 
cored halo profiles evident in many galaxies, and may yield 
more realistic substructure than collisionless cosmogonies. 
Early numerical studies of weakly self-interacting halos 
mimicked SIDM using particle codes with Monte Carlo scat- 
tering, which bred realistic cores but found a gravothermal 
catastrophe that could eventually degrade the cores in to 
isothermal cu s ps (IMoore et al.ll2000l; lYoshida et al]|2000allbl ; 



isothermal cusps (Moore et al. 2000; Yos 
lBurkert|[2000l ; lKochanek fc White! |2000| 



... iDave et alj|200ll) 
Howe ver later analyses ( Balberg et al]|2002i ; Ahn fc Shapiro! 



120051 ) considered more general initial conditions, delaying 
collapse beyond the Hubble time. Polytropic halos (as in 
this work) can describe a more strongly interacting fluid 
regime, where the mean free path is short enough that 
conduction and gravothermal effects vanish. 

Despite the indications from galaxy scales, the fash- 
ionable preference is to defer and displace the faults of 
CDM substructure onto "baryon feedback", which is be- 
set with long-term challenges in numerical methods and 
in theory. The evidence on cluster scales is still ambigu- 
ous enough to allow this. Many observers assume cuspy 
profiles in their data fits. (Fully non-parametric modelling 
remains rare.) For relaxed clusters, some X-ray deprojec- 
tions show soft-cored halo profiles, (e.g. | Nevalainen et al] 



19991; lEttori et alj l2002bl; iKatavama fc Havashidal |2004| ; 
Voigt fc Fabianll2006l ; IZhang et alj 2005, 200lj) while others 



seem compatible with cusps (e.g. Pointecouteau et al.ll2005l : 
IVikhlinin et~aT]|2006l ). In some cases the total mass profile 
appears cuspy, but not to the ex tent expected of a coll i- 
sionless halo affected by gas (e.g IZappacosta et al] I2OO6I ). 
Gravitational lensing analy ses also give mixed signs: some 
prefe r or allow soft cores (|Tvson et ail Il998l; iDahle et al " 



Gayazzi et all 120031; iDiego et alj 120051; ISand et al 



2002 . |2004 120081 ; lHalkola et aljbood |200S|; iRzepecki et al 
2007) while others prefer cusps (IBroadhurst et alj |2005[ 



Sharon et aT]l2005l ; ISaha et alj|2006l ; iLimousin et alj |2007) 



In more violent circumstances, cluster mergers have 
been treated as probes of dark interactivity. One gravita- 
tional lens "bullet cluster" was claimed as a merger of col- 
lisio nlessly interpenetr ating halos, separating from shocked 
gas (|Clowe et alj|2006l ). The mass, speed, timing and rarity 
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of the hypothetical merger have be en questioned, and well- 
tuned simulations dev i sed in reply (lHavashi fc White] | 2006 |; 
Farrar fc Rosenl 120071: Izhad 120071: ISpringel fc Farrarl I2OO7I; 
Milosavlievic et all |2007| ; lAngus fc McGaughl |200S| ; iNusserl 
20081 ). Subsequently, an anti-bullet cluster has appeared, 



where the dark matter is the more dissipative constituent 
(settled in the middl e) while galaxies fly on the periphery 
jMahdavi et all 120071 ) . Another lensing cluster is reported 
with an encircling ring or shell of dark matter (|jee et al.l 
120071 ). Taken together, these special cases tell an incon- 
sistent story about dark physics. However, particular pro- 
jected morphologies admit multiple interpretations: for in- 
stance the "bullet" subcluster velocity vector can be re- 
versed, and the line-of-sight shapes and displacements are 
unknown. This paper cannot aim to disentangle all the latent 
assumptions in the dark matter merger problem, but clearly 
some alternative gestalts are needed. Varieties of SIDM re- 
main among the promising candidates. 

In the absence of a central mass or gas inflow, a poly- 
tropic halo can have a constant-density core (Appendix lD2[) . 
compatible with galaxian evidence. However for galaxy clus- 
ters with inflowing gas, we do not obtain simple cored pro- 
files like those assumed commonly. The presence of inflow 
requires a central mass for stationarity, and a polytropic 
cluster halo grows a density spike within its core. The spiky 
halos are effectively a "contraction" induced by the central 
mass and gas inflow. For large F2 or large m, the core ra- 
dius may be small enough to give a misleading appearance of 
a N FW-like cusp. Sma l lness of some observed cluster cores 
(e.g. lDahle et al.ll2003l ; iKatgert et al 1 12004 iLimousin et al] 
120071 ) is not evidence against SIDM. Rather, it will help 
to constrain the dark freedom F2 and the inflow history of 
baryons. Cluster cores as small as a few tens of kiloparsecs 
are possible if F2 is large. This fits the concordance favour- 
ing 8 < -F2 < 10, which minimises the cent ral mass m* (this 
work) and agrees with iNunez et al] (|200rj ). who effectively 
find F2 ~ 9.6 by fitting galaxy rotation curves. 

Given any alternative closed set of equations for the 
dark dynamics and statics, one can repeat the formulation 
of this work, to find obtain another set of differential equa- 
tions coupling the gas and halo. The interplay of these con- 
stituents in their shared potential must always lead to the 
exclusion of some domains due to cooling catastrophe or 
acoustic breaks in the gas. However it must be proven, for 
each scenario, whether nonzero m* is required (as in the 
model dipolytropes here) and whether structures with cos- 
mic baryon fraction can exist. 



5 CONCLUSIONS 

We have presented a self-consistent 2-component model for 
galaxy clusters, bound by a non-static gravitational poten- 
tial that emerges naturally from the solutions along with 
the co-adapted halo and gas profiles. Applying this formu- 
lation to clusters of plausible total mass and composition, 
we reconcile some of the observational difficulties involving 
gas inflows due to cooling. Furthermore, we find that sta- 
tionary solutions of the cluster structure invariably require 
(or develop) a non-zero central mass. 

We have analysed the distribution of cooling gas in a 
responsive spherical halo of megaparsec scale. Mass, mo- 



mentum and energy continuity arc imposed. Bremsstrahlung 
radiative cooling is allowed to become dynamically signifi- 
cant. All constituents participate gravitationally. Realistic 
models emerge when dark matter has a polytropic equa- 
tion of state, which is justified in terms of the equilibria of 
Tsallis thermostatistics, adiabatic dark self-interactions or 
Bose-Einstein condensation. 

We find that there exist steady, continuous solutions 
spanning all radii inside the halo. The joint constraints of 
cooling and acoustic continuity set the minimal central mass. 
The minimal m* varies with m and F2 but only weakly 
with the gas surface temperature Tr. The cluster's total 
gas fraction is linked with rh, R and Tr. The masses, rh, 
densities, temperatures, velocities of any particular solution 
can rescale to yield another valid solution with the same 
radial dimensions. 

A cosmic baryon fraction and observed SMBH masses 
are consistent with the halo's effective microscopic degrees of 
freedom being in the interval 7 < F2 < 10. The lower limits 
on m* are laxer if m 2 /m 3 is smaller. For cosmic composition, 
the fiducual total cluster mass and rh ^ 1 m© yr _1 , we 
always find m, > 5x 10 s m Q (or m*/m > 2 x 10~ 9 ). Smaller 
central masses are impossible in steady clusters, unless extra 
physics dominate. To enable m* as small as 10 6 -10 7 m@ m 
a cluster, we need 9 < F2 < 1 . This agrees with galaxy 
rotation models of I Nunez et al.l l|2006h which imply F2 « 
9.6. 

The halo density develops a spike around the central 
mass, surrounded by a flat core attenuating to a fringe on 
megaparsec scales. These layers are less distinct when F2 
is larger. This varied structure remains apparent in projec- 
tions of the total column density. For F2 high enough to 
enable plausible m„, we find halo cores with 10 kpc < Ri < 
300 kpc. Observationally, th ere are reported core fits with 
15 kpc < R x < 200 kpc (e g. iDahle et al.|[200^ ; iDiego et al] 
l2005l ; lRzepecki et alj|2007l ; lrlalkola et al]|2008l ) which is sim- 
ilar to what we obtain. We predict that mass profiles steepen 
beyond NFW in the dim fringe. Given observable scale radii 
such as R2 , one could predict the outer radius where a halo 
naturally truncates. 

Our solutions belong outside the inevitably overcooling 
regime where classic cooling flows were constructed. They 
naturally provide a non-zero floor temperature, obviating 
the need for (unobserved) mass dropout and cold condensa- 
tion. The entropy, density and temperature profiles broadly 
resemble observed clusters, suggesting that varying the gas 
parameters may enable detailed fits in future. We find a 
shallow gas entropy ramp at radii inside the sonic point, 
rather than a flat isentropic core. In projection, the intra- 
cluster medium resembles the classic /3-model X-ray profile 
in the outskirts, plus a peaked cooling/warming core. The 
central luminous spot is comparable to AGN power, though 
opacity, conduction and detailed AGN anatomy could prob- 
ably soften and spread the emission, and lower the radiative 
efficiency considerably. In the optimal-m» solutions, the ra- 
tio of peak/dip temperatures is a factor 5-40. This reduces 
the need for AGN self-regulation (as distinct from incidental 
heating). The inclusion of thermal conduction and Fi > 3 
might improve the profiles and lessen the AGN role further. 

By construction, the classes of solutions that 
we obtained comprise the steady configurations of 
bremsstrahlung-cooling Mpc-scale spheres. If our broad 
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physical assumptions hold, and if structural asymmetries 
are subdominant, then these solutions represent end-points 
of cluster relaxation. In the configuration-space of spher- 
ical clusters, the solutions are fixed points. If feedback 
or non-gravitational heating processes become globally, 
persistently important, then they drive the evolving cluster 
state in a forced orbit around those natural points. 

A temporary disturbance of the system may cause a 
local cooling catastrophe or acoustic disconnection, and 
the structure must somehow adjust until reaching another 
steady state in neighbouring configuration-space. Our anal- 
ysis does not outline a particular evolutionary path, but the 
profiles offer clues. The critical bottlenecks for gas continu- 
ity occur at radial scales typical of elliptical galaxies. This 
may be a preferred layer for cold gas dropout and star for- 
mation, during any transient, externally driven detour into 
cooling catastrophe. Stability analyses are needed to deter- 
mine whether our scenario implies regulation or runaway 
monolithic collapse. For F2 > 6, the innermost halo is only 
marginally Jeans-stable, implying that large-amplitude dis- 
turbances could trigger a local gravitational collapse of dark 
matter (without involving the gas directly). This mechanism 
for dark growth of SMBH may turn out to be an important 
process (besides baryonic feeding and gravitational ejection) 
influencing SMBH demographics. 
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APPENDIX A: NATURAL UNITS 

If the gas has approximately solar composition then 
the bremsstrahlung constant has a value B = Bq « 



5.06 x 10 1 



iRvbicki & Lightma 


n 


|l979|) 


as in Saxton et al 




bood) 


lAnders & Grevesse 




1989h. 



calculated according to 
r a Gaunt factor gB = 1.25) 
with the abundance tables of 
We parameterise the composi- 
tion dependency of B relative to the solar value as a cor- 
rection factor, ( = B/Bq. Values of £ depend on the 
abundance- weighted ionic mass (m), charge (Z) and charge- 
squared (Z 2 ), 



3/2 



(Al) 



Z \m/m c + Z 

For astrophysical plasmas, £ is close to unity: for the same 
<7b in a pure H plasma, £ ~ 0.979; for a 9 : 1 mix of H and 
He, C « 0.925. 

If we define a system of units such that B = 1 and the 
gravitational constant G = 6.6732 x 10~ 8 g _1 cm 3 s~ 2 = 1 
also, then the unit of length is 



U x 



B/G = 7.58 x 10 23 C cm = 0.246 ( Mpc , 



(A2) 



It may be significant that this scale, which is natural 
to any object governed by self-gravity and optically-thin 
bremsstrahlung, is typ ical of the observed X- ray core radii 
of galaxy clusters (e.g. Ijones fc FormarJ[i"984h . 

If we adopt a velocity scale where a 2 — 1 corresponds 
to a temperature of 1 keV, then the unit of velocity is 

U v = 3.95 x 10 7 cm s" 1 , (A3) 

which implies a time unit 

U t = U x /U v = 1.92 x 10 16 C s = 0.608 C Gyr . (A4) 

The age of the universe (|Spergel et all l2003h is presently 
thought to be ~ 22.5/C- The unit of mass arises from 

1.77 x 10 46 C g = 8.91 x 10 12 C m e , (A5) 



and the unit of density is 

U p = U m U~ 3 = 4.07 x 10" 26 C" 



g cm 



(A6) 



The critical density for the universe today is p c = 2.33 x 
10~ 4 C 2 U p , and the mean matter density is fi m p c ~ 6.23 x 
10~° £ 2 Up. The unit of particle number density depends on 
plasma composition in a more complicated way, 



Ur, 



U 



m + Zm 



(A7) 



For solar composition we have U n w 3.96 x 10 2 cm 3 . The 
units of mass accretion and power are 

U m = U m /U t = 1.47 x 1O 4 M yr" 1 , (A8) 
U L = U m UllU t = 1.44 x 10 45 erg s" 1 = 3.77 x 10 11 L .(A9) 
and neither depends on composition, £. 

APPENDIX B: SCALING RELATIONS 

Given one steady cluster model, it is possible to construct 
a set of equivalent models that differ only by uniform mul- 
tiplicative rescaling of the physical variables. Let us define 
the transformation factors as 

r — > X r r, 
m — > X m m, 
m — > Xm m, 

2 v 2 



v — > X v v and 
X p p. 

Mass conservation © implies a constraint 



Xr, 



X r X p X v 



(Bl) 
(B2) 



The equations of the mass profile, such as (|15|) or (|4ip . re- 
quire that 

Xm — XpXr . (B3) 

Mach numbers must be left unchanged during the transfor- 
mation, and thus 



X,, — X a 



(B4) 



In each of the gas equations, say (|31|l . all of the additive 
terms must rescale by the same product. After some evalu- 
ation, this implies that 

> 



X 



(B5) 



Satisfaction of (|B3|) and (|B5|) implies that the spatial dimen- 
sions cannot vary, 



X T 



(B6) 



Therefore, any valid similarity transformation parame- 
terised by a scale X, implies the following scaling factors 
for the key physical variables: 

Xm = X p = X a = X , 

x v = x 1 / 2 , 

J^ m — -A. , 

X s = X (F ~ 2)/F . (B7) 
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Luminosities and surface brightnesses scale as Xl = 
XrhXy = X 5 / 2 , so the X-ray luminosity scales as m 5//2 
within any family of solutions. The relations (|B7[) imply the 
existence of two invariant length scales, associated with the 
mass inflow and temperatures: 



Rri 



Ra 



7K7im/m|j 



(B8) 
(B9) 



The latter is algebraically equivalent to the sonic radius in 
simple, adiabatic Bondi accretion, although cooling and self- 
gravity mean that our models needn't develop a sonic point 
nor M 2 extremum at this radius. Together, the parame- 
ters (F\,F2,R, Rm, Ra) uniquely denote a set of homologous 
cluster models. 



APPENDIX C: 
FREEDOM 



EFFECTIVE DEGREES OF 



The key property of an ideal fluid is F, the effecitve degrees 
of freedom. Thus F is a qualitatively decisive parameter 
of the halo models. In ordinary space, free particles have 
three translational degrees of freedom. If the particles are 
individual, point-like entities lacking substructure then F = 
3. However many physically motivated scenarios entail F > 
3 or non-integer values. 

If individual particles can rotate, twist or distort then 
there are additional internal microscopic degrees of free- 
dom (e.g. F = 5 for a diatomic gas). Highly relativistic 
or radiation-dominated fluids have F = 6 (e.g. cosmic ray 
contributions to ICM pressure). Larger integer values of F 
could also occur if the particles experience higher spatial di- 
mensions, e.g. if their de Broglie wavelength is smaller than 
the scale of compact hidden dimensions. If a fluid includes 
subspecies that do not fully interact, then the effective F is 
larger than for single species. 

Some alternative scenarios involve fewer degrees of free- 
dom. If the fluid is a classical Bose-Einsten condensate (as in 
Goodman! I2OO0I ; lArbev et aill2003l ; iBohmer fc Harkoll2007l ; 

A 



2008) then F = 2, for an equation of state p oc p 2 . 



case F = 1 could describe constrained particles, analogous 
to beads on an abacus. An incompressible fluid corresponds 
to F = 0. Isobaric conditions can be described by F = —2. 

If non-local physical interactions are important, a nd the 
medium is described by Tsallis' statistics (Tsallis 1988), then 
F is effectively so me non-integer, F = (3q — l)/(q — 1) for 



some constant a JPlastino fc Plastino 1 19931 ; lHanse r\ 120051 ; 
iNunez et alj [20061 ; IZavala et al.l 12006 ). In a gravitational 
context, this includes and entails the ephemeral constraints 
and interactions present on all intermediate levels between 
the small scale of two-body scattering and the large scale of 
the global potential. Other meso-scale physics, such as the 
energy associated with turbulent eddies, can also provide 
larger, non-integer values of F. If highly efficient heat trans- 
port processes operate then the fluid approaches isothermal- 
ity and F —* 00. 



APPENDIX D: STANDARD GASLESS HALO 
MODELS 

Dl comparative measures of a halo 

To found our treatment of two-fluid cluster models, we will 
here review the intrinsic properties of polytropic halos (with- 
out gas), and contrast them with other halo models in the lit- 
erature. In order to compare theoretical and semi-empirical 
halo models with each other and with observations, it is nec- 
essary to define some global physical measurements. Every 
real halo ought to have a finite outer radius, R, but the 
invisibility of dark matter means that R isn't directly deter- 
minable. In practice, the cluster's baryons are only visible 
out to certain detection limits (e.g. to a limiting X-ray flux) 
and this extent sets lower bounds on R. We prefer to charac- 
terise the models using spatial measurements, (which can be 
compared to the true outer radius) since these are invariant 
under mass rescaling (Appendix [Bj. 

For non-singular and cored halo models, a "King ra- 
dius" is defined in term s of the central conditions, (e.g. 
iBinnev fc Tremaindfl987l . p. 228) 



Rk = ^9a 2 /ATvGp\ 



(Dl) 



This scale typifies the extent of the flat density core in many 
models. The cluster rescaling (appendix[B} leaves _Rk invari- 
ant. However Rk is undefined for cuspy halos or those with 
a central point mass. Thus we require alternative measures 
of core size and overall halo concentration. 

Firstly, let us define a radius that contains a majority 
of the mass, or that typifies the central concentration. We 
refer to a sphere's total mass m and moment of inertia, 



8tt 



pr 4 dr 



which give a mass-weighted lever radius, 



Ri = s/M/2m 



(D2) 



(D3) 



which is scaled such that Ri — R for a uniform sphere. 
The radius Ri is applicable to models where the core is 
not explicitly parameterised. Since density decreases mono- 
tonically in r, the inner layers dominate Ri. When the the 
mass is centrally peaked or the core is small compared to 
the fringe, the ratio Ri/R is small. Though the true sur- 
face may be invisible below some flux or density thresholds, 
truncated observational estimates of Ri or Ri / R might still 
approximate the global values acceptably. 

We define another radial scale measuring the cluster's 
self-gravity. The gravitational potential energy of a spheroid, 



W 



-4ttG 



mpr dr = 2tv / $pr dr 



f 

Jo 



(D4) 



is finite for realistic models. This leads to the definition of a 
gravitational radius, 

R w = ~Gm 2 /W , (D5) 

fsee lBinnev fc Tremainelll987l . pp. 68-69.), which can be re- 
garded as the size of energetically typical orbits in the halo. 

Halo structure can be diagnosed by observable kine- 
matic tracers, such as velocities of gas clouds or stars in 
circular orbits in a galaxy, or the motions of galaxies within 
a cluster. The circular orbital velocity of test particles peaks 
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at some radius R if 



(Airpr -m)| r=J?o 
d In p 



= , 

< -2 



and 



(D6) 



We denote a sequence of signature radii where the log- 
arithmic density slope passes specific values: 



din p 
dlnr 



r={K 1 ,H 2 ,«3,«4,-} 



= {-1,-2,-3, -4,...} 



(D7) 



For instance, p oc r~ 3 at the slope-3 radius Rj,, and so on. 
These slope radii may be multi-valued, if the density profile 
undulates (i.e. exhibiting concentric, alternating steep and 
flat layers). If the halo is radially finite (R < 00) then all 
slope radii are finite too. Infinite models may lack some of 
the slope-radii. If the rotation curve peaks anywhere, then 
(|D6j) implies that R ^ R2. The total mass cannot be finite 
without 7?3 existing. A finite moment of inertia and Ri re- 
quires finite i?5. Finite W and R m require R2.5 and -R3 to 
occur at least once in the outskirts. 

The literature on cosmological simulations convention- 
ally defines a "virial radius", R v , enclosing a mean density 
that is some multiple, 5 C , of the cosmic critical density. This 
overdensity is in the range 10 < S c < 200, dependin g on 
the cosmological model. As in iBrvan fc Normanl dl998t ). we 
use 5c = 187T 2 + 82(f2 m — 1) — 39(f2 m — l) 2 , corresponding 
to idealised spherical collapse, and f2 m = 0.27 such that 



97.01p c « 0.02259 U p 



(D8) 



3m v 
4-kR* 

The virial radius conveniently measures idealised, radially 
infinite models, or numerical simulacra which are unresolved 
in their fringes. The drawback of R v is that it loses infor- 
mation about the outskirts. It also fails to characterise com- 
pact objects where R v encloses the entire mass (e.g. the 
most compact (ra*, T)-optimal models with R < 1 Mpc). 
Since R v is defined relative to an absolute density, it does 
not transform neatly under mass and temperature rescaling 
(Appendix [B]). On the other hand, it is always possible to 
rescale the cluster masses so that R v = R or some other 
signature radius ( H3.5|l . The same is true for King models or 
any model with at least one free scale. 

X-ray imagery and gravitational lensing studies con- 
strain the column densities of gas and dark matter, as pro- 
jected onto the plane of the sky. Thus it is useful to calcu- 
late comparative two-dimensional projected properties. An 
effective half-light radius, R c , is conventionally defined by a 
line-of-sight cylinder that encircles half the emission (or pro- 
jected mass). If the total mass and halo radius were known, 
then we can also define a radius ife for the image contour 
with mean brightness (or column density), £ = m/ivR 2 . 



D2 gasless polytropic halo 

Here, for the sake of clarifying our main results, we review 
the intrinsic behaviour of finite polytropic dark halos without 
gas. We will examine the influence of a central point-mass. 
The upper blocks of Table ID1I characterise a set of gasless 
polytropic halos, of different degrees of freedom, with and 
without a central mass. 

The first subset lack a central mass (ra, = 0), and share 



the same entropy and central density: s = 1, p(0) = 1. 
Figure ID 1 1 shows their density profiles. These s olutions are 
classical Lane-Emden spheres of index n = F/2 (|Lanelll8"7(il ; 
lEmden|[l907l : IChandrasekharll 19391 : iHoredtll 198^ 1. We prefer 
F as the more physically motivated notation. Each sphere 
has a core of nearly uniform density, surrounded by declining 
outskirts. For smaller F the core is a larger fraction of the 
volume, and the fringe is steeper. If F < 10 then the halo 
possesses a zero-density outer surface at radius R. These 
finite polytropes don't tend to any asymptotic outer slope; 
the density index steepens infinitely as r — > R. 

With central conditions held constant, the radius and 
mass increase with F, and the binding energy increases both 
in absolute terms and per unit mass. In absolute terms, Ri 
increases slightly with F. Proportionally, the core (Ri/R), 
the gravitational radius (R w /R), the two-dimensional radii 
(Rs/R, R c /R) and the rotation peak (R a /R) all shrink with 
F. The slope-radii (-R2, R3, Ri) also shrink with rising F. 
They become multi-valued for large F (the index of p wob- 
bles in some layers). In such cases, we tabulate the steep- 
ening point farthest on the edge of the core. This tends to 
be near the peak of the rotation curve. We generally find 
that R > R3 and Ri > R4. Consecutive slope indices 
(R2, R3, Ra ■ ■ •) occur in roughly even steps, but the steps ul- 
timately diminish near the true surface (R) . In non-singular 
halos, J?™ ~ Ri. As F increases, all the signature radii shrink 
relative to R in a common manner (see Figure |D2)| . 

However the non-singular model is a specially contrived 
condition. Most galaxies are thought to contain a point-like 
central mass such as a nuclear star cluster or black hole. 
At cluster scales, the analogous object is a cD galaxy or its 
black hole. To clarify the effect of such a mass, we tabulate 
singular polytropic halos with central masses in astronom- 
ically realistic proportions: m«/m = 10 -8 , 10 -6 , 10 -5 . We 
fixed the total mass and R of corresponding Lane-Emden 
spheres, but vary the entropy. Figure |P3l shows density pro- 
files for m»/ra = 10 -8 . The central object draws a density 
spike about itself, of index —F/2. Beyond this sharp sphere 
of influence, the halo flattens into a core, then steepens into 
outskirts and a surface like those of nonsingular models. 
For small or medium F, the addition of ra* > reduces 
the slope-radii (R2, R3, Ri) slightly. However for sufficiently 
large F and m„, the spike steepens the entire core, to the 
point where R2 and higher slope-radii vanish. 

For F > 9, the spike develops density undulations (at 
sub-parsec scales for a cluster). In this tiny, deep core, the 
slope radii are multi-valued and the rotation curve peaks 
multiply. Some undulations locally approach the brink of 
Jeans stability. The dense spike reduces Ri and R w dramat- 
ically for F > 8. 

Table [DTI reveals several trends for as m*/m varies. As 
ra, increases, R w /R changes appreciably before the slope- 
radii and rotation peak do. The two-dimensional projected 
quantities (Rs/R, R c /R are the least sensitive to the central 
mass. Large F enhances the sensitivity of all the of signature 
radii with respect to ra*. 

If F ^ 10 then the halo density attenuates indefinitely 
(R = 00), regardless of ra». The borderline case of F = 10 
has an infinite r adius but finite mass: this is the well-known 
iPlummerl (|l91ll ) model, where p oc r -5 at large radii. In the 
isothermal limit, F — > 00, the fringe declines like p oc r~ 2 
( Appendix ID3|) . 
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Figure Dl. Density profiles of polytropic dark halos with dif- 
ferent adiabatic indices, 7 = 1 + 2/F, but no gas and no central 
mass. We set s = 1 and p = 1 at the origin. Labels denote the ef- 
fective degrees of freedom F in each case. For large F the fiat core 
fills less of the total volume. The case F = 00 is the non-singular 
isothermal sphere ("NIS", qD3|) . 




4 6 
degrees of freedom F 

Figure D2. Signature radii (as fractions of the surface radius 
R) of non-singular (m* = 0) finite polytropic halo models. Bold 
curves mark the density slope radii, R2/R, R3/R, R4/R and 
R5/R. "I" marks the lever radius, Rj/R. "M" marks the half- 
mass radius, R m /R. "W" marks the gravitational radius, R w /R. 
"O" traces peaks of the rotation velocity, R / R. Two-dimensional 
projected quantities are marked "£" (R^/R the mean- brightness 
radius) and "e" (R e /R the half-light radius). These radii stay 
bundled together but shrink in relation to the surface as F rises. 



D3 Isothermal spheres 

If the particle velocities are isotropic and a 2 is constant 
everywhere, then the halo is an "isothermal sphere." This 
is essentially an extreme polytrope in the limit F — > 00. 
Isothermality is plausible when some strong mechanism as- 
serts global thermal equilibriu m: e.g. thorough and violent 
relaxation (|Lvnden-B cll 1967), or efficient thermal conduc- 
tion. Isothermal spheres are popular toy models in gravita- 
tional lensing studies. Such halos also predict flat rotation 
curves, resembling the observed outer parts of disc galaxies. 

The density profile depends on the central boundary 
condition. The "non-singular isothermal sphere" (NIS) has 
dp/dr = at the origin, exhibiting a shallow density core, 
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Figure D3. Density profiles of gasless polytropic halo models in 
the presence of a central gravitating mass (m* = 10~ 8 m). Curves 
are annotated with their respective F values. 



but at large radii i t ten ds to a decline oc r 2 (see e.g. 
iBinnev fc Tremainel fl987l . and the F = 00 curve in Fig- 
ure [Dl}. The "singular isothermal sphere" (SIS) has a self- 
similar profile, 



p = a 2 /2-KGr 2 



(D9) 



For NIS and SIS, R3 and higher slope radii never occur. The 
density fails to vanish at any finite radius, so the halo lacks a 
distinct surface. Within some ad hoc truncation radius, the 
mass, moment of inertia and gravitational potential energy 
of the SIS are 



m = 2ar/G , 
I = 4aV/9G 

W = -4cr 4 r/G 



and 



(D10) 



It follows that the gravitational radius R w — r and the ef- 
fective lever radius Ri = y/Er/3. For the SIS, the circular 
velocity is radially constant. A SIS gravitational lens bends 
light rays by a constant angle at all projected radii. The SIS 
virial radius is 



R\ 



3a 2 



1/2 



4.598a 



(Dll) 



2%G8 c p, 

which is w 1.130 Mpc for a halo at 1 keV temperature. 
D4 pseudo-isothermal sphere 

In some observational studies the exact non-singular isother- 
mal sphere is approximated by an empirical cored pro- 
file that also has oc r~ 2 outskirts. The "pseudo-isothermal 
sphere" (PIS) has a density profile 

where x — r/r s and (p s ,r s ) are some density and radial 
scales. The PIS has infinite radius, mass and moment of 
inertia. At large radii, Ri/r — v5/3 and R m /r — 1. The 
slope-2 radius is at infinity, as is the peak circular velocity 
(R2 = Ro = 00). Figure [D4l shows numerical solutions for 
the virial radius depending on (p s ,r s ), comparing PIS to 
some other radially infinite models. 
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D5 Hubble profile and /9-model atmospheres 

The modifield Hubble profile (|Hubbld Il930l ; iRood et al.1 
ll972l : lKindll972h approximates the projected mass or bright- 
ness of any non-singular cored sphere as 

E ~irb' ( D13 > 

defined such that the intensity drops to half its central value 
(So) a t the projected radius scale r B . This empirical fit has 
been applied to the stellar light profiles of elliptical galax- 
ies, globular clusters and galaxy clusters. The corresponding 
spatial density profile is 

This distribution yields more realistic rotation curves than 
the isothermal mode ls, as there is a peak at fi nite radius, 
R « 2.920r s (e.g. iBinnev fc Tremaine] Il987l . pp. 39-41). 
However the model is only applicable in and near the core, 
otherwise it implies infinite radius, mass, potential energy 
and moment of inertia. As r — > oo we have Ri — > oo, 
R w — > oo, Ri/r — > and R w /r — > 0. The density slope- 
radii are R2 = y/2r B and R3 = 00. 

The model was extended to describe the X-ray emitting 
gas by assumi ng isothermality and local hydrostasis in clus- 
ter cores, e.g. iLea et all (|l973T ) who assigned identical tem- 
peratu res to gas and galaxies. ICavaliere fc Fusco-Femiano! 
ljl976f l assumed a galaxy/gas temperature ratio of (3 — 
crl/ at, deriving a gas density profile 



10 1 F 



P = pa [1 + (r/r c ) 



-3/3/2 



and X-ray surface brightness profile 

s = so [1 + (b/b c ) 2 r 3P+i 



(D15) 



(D16) 



where b is the projected radius. This became a com- 
monplace fitting forumula for X-ray imaging observa- 
tions, (e.g. Bahcall fc Sarazin 19771 ; Gorenstein et al] 1978 



Branduardi-Ravmont et al.l Il98ll: IJones fc Formanl 1984 



Neumann fc Arnaudll 19991 ). Shorn of the original isothermal 
assumption, more recent studies adopt either f|D 15f) or ()D16[) 
as a conventional template, and seek to infer actual radial 
variations of temperature from the data. It is worth empha- 
sising that the /3-model is a parameterisation or an ideal- 
isation of the innermost observable regions. The ultimate 
outer density index is —3/3, and the implied asymptotic gas 
mass is infinite unless /3 > 1 (according to (|D15[) 1. However 
()D16|) implies that the luminosity is infinite if 1 < /3 < 2. 
This breakdown implies that the density and/or tempera- 
ture must attenudate even more steeply in the outskirts of 
real clu sters. Indeed, some observations show that f3 steep- 
ens (e.g lVikhlinin et~al ] |l999l : lNeumann|[2005h . 

D6 Hernquist profile 

iHernauistl (|l990i ) proposed an analytic model for galaxy 
spheroids, with density profile that attenuates infinitely 

(D17) 



27nr(l + x) 3 r| 
but nonetheless yielding a finite mass profile 
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Figure D4. The dependency of the virial radius upon the halo's 
density at its scale radius, for: pseudo-isothermal NFW (N); 
Burkert (B); Hernquist (H); Sersic n = 2, 4 (S2,S4) models. 



The moment of inertia is infinite. The density slope passes 
integer values at X2 — |, x% — 2 and xa = 00. The rota- 
tion curve peaks at the scale radius, x a = 1. The projected 
central brightness is infinite, but the outskirts decline fast 
enough that R c is finite. Figure ID4I shows the relation be- 
tween virial radius and (p B ,r s ). 



D7 N-body simulacra 

iDubinski fc Carlberd (|l99lh found that collisionless struc- 
tures emerging in N-body cosmological simulations develop 
a cus py, power-law central den sity profile. The "NFW pro- 
file" (|Navarro et al.|[l996l , Il997h is a popular empirical fit to 
such halos, with a density that follows 



P = 



x(l + x) 2 



(D19) 



where x = r/r s and (p s ,r s ) are fitting parameters of a 
particular halo. These parameters follow trends in rela- 



tion to the halo m a ss, which depend o n cosmology (e.g . 
Bullock et all l200ll : iDolag et al.1 12004 IShaw et all l200d : 



Maccio et al.|[2007l ). Physically, the scale radius r B = R2, 
the slope-2 radius. This radius defines a concentration, 
c = R v /r s . A halo truncated at some finite radius has a 
mass and moment of inertia of: 



m = 4irp s r B 



8tt 



-p*r B 



ln(l + x) 



3 1 
2 + 1 +x 



(D20) 



+ 



(a:+l)(x-5) 



+ 31n(l + a;) 



.(D21) 



There isn't an outer surface, and both m and I are infinite as 
x — > 00. The inertial radius Ri — + cxj at infinity, while the its 
concentration ratio vanishes (Ri/r — > 0), which means that 
the rotational properties depend on an ad hoc truncation 
radius. The gravitational potential energy is finite, 



W 



-8n 2 Gp B r B 



(D22) 
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Thus R w — > oo and R w /r — > at large radii. The rotation 
curve peaks at R w 2.163 r s . Using ([D19)) and (fD20l) . the 
virial radius equation (|D8|l is transcendental. (See Figure lD4l 
for numerical solutions. ) 

iMoore et all <|l99Sh and IZhad (|l996h proposed varia- 
tions and generalisations to the NFW formula, consisting 
of a broken radial power-law again, but with different in- 
dices. Expressions for the global quantities differ slightly 
from those above, but the models are qualitatively similar: 
infinite in mass and radi us, and ill-defined rotational prop- 
erties . More recent work (jMerritt et al.|[2005l ; iGraham et al.l 
I2006T ) suggests that Sersic profiles (Appendix ID lOf) fit sim- 
ulated collisionless halos better. 

Density slopes of the simulacra are least certain on the 
outskirts and near the origin, both places where mass res- 
olution degrades. Cuspy profiles emerge consistently from 
cosmological N-body simulations, but the causes of this 
shape still lack a comprehensive analytic derivation. Possi- 
ble causes may involve: cosmic expansion, ongoing infall and 
accretion history; the simplifying assumption of a collision- 
less medium; the approximation of discretised mass; implicit 
low-pass filtration in numerical Poisson solvers; or perhaps 
other less obvious computational and physical factors. 

D8 Burkert profile 

Burkcrt ( 1995) proposed an empirical halo model based on 
observed rotation curves of halo-dominated galaxies, 

Ps 



0.6 



P = 



(l + x)(l + x 2 



(D23) 



where we denote a normalised radius again, x = r/r B . The 
density index is —2 at R2 fa 1.521 r s . The index ultimately 
approaches —3; R3 — 00. There is a flat density core, like 
in the nonsingular polytropic and King models (and unlike 
the cuspy NFW and Sersic profiles). The mass and inertial 
moment enclosed at x are 



m = Trps^s [2 ln(l + x) + In (l + x 2 ) — 2 arctanir] 



2tt 5 



2x — 4x + 2 arctan x + In 
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and both are infinite as r — > 00. As with NFW, the Burkert 
halo mass is intensely centrally concentrated: Ri/r — > as 
r — * 00. The gravitational potential energy is finite, 

W -- 
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and R w /r — > at large radii. The rotation curve peaks at 
R ~ 3.245r s . Like the NFW halo, the virial radius equation 
is transcendental (Figure [D4] shows numerical solutions). 

D9 King model 

King (1966) presented a cored stellar-dynamical model, de- 
rived from first principles. Its basis is a phase-space density 
function, 



-a(* + u 2 /2) 



-a*t j 
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assuming locally isotropic particle velocities, truncating at 
some escape energy corresponding to the equipotential ($ t ) 
of a zero-density outer surface at "tidal radius", R = r t . 
The model is a self-consistent description of a non-isolated, 
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Figure D5. Signature radii relative to the tidal surface, for King 
models with various concentrations, c = log 10 (ii/ Rk)- Lines and 
annotations are the same as in Figure ID2I The signature radii 
span much of the halo volume (unlike high-F polytropic halos). 
Their ratios stay roughly steady as c increases. 



self-bound, collisionless sphere. It wa s originally app l ied to 
globular clusters with escaping stars. iFirmani et al.l (1200 lh 
reapplied it to cored, thermal, self-interacting dark halos. 
The mass, moment of inertia and gravitational potential en- 
ergy are all finite. The local density, 



depends on the dimensionless potential offset, 



; [$ t - $(r 



(D28) 
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and F is the lower incomplete gamma function. At the outer 
boundary, p = ip = 0. Radial coordinates for the equipoten- 
tials are obtained by solving the Poisson equation, subject 
to the inner boundary conditions ip > and Vtp = 0. The 
latter condition precludes a central mass, m* = 0. 

Like the F < 10 polytropes, the King model is ra- 
dially finite. Thus the slope radii (R2, R3, R&) are also fi- 
nite, and the density index attains large negative values 
near the edge. Concentrations are conventionally denoted by 
c = log 10 (i?/i?K). Figure |P5l shows the variation of signa- 
ture radii with c. Both R2/R and R /R drop rapidly with in- 
creasing c, and they are multi- valued for highly concentrated 
models, c > 2.7. The other main signature radii vary only 
within factors of a few in the shown domain (0 ^ c ^ 4.5). 

D10 Sersic profile 

The ISersid (| 19681 ) model is an empirical fit to the two- 
dimensional projected starlight of spheroids such as elliptical 
galaxies and spiral bulges. 



E = E e 
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The radial coordinate is scaled in terms of the half-light 
radius, x = r/R c . The shape parameter n ~ 4 for elliptical 
galaxies (the classic profile of lde Vaucouleurslll948r i or n ~ 2 
for galaxy clusters. The parameter b depends on n implicitly, 
via lower-incomplete and complete gamma functions, 



2r(2n, b) = T(2n) 



(D31) 
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ICiotti fc Bertinl (|l999h derived a series expansion, b ~ 
2n - 1/3 + 4/405n + 46/25515n 2 + 131/1 148175n 3 - 
2194697/30690717750n 4 . Sersic profiles appear ubiquitous 
among stellar spheroids in nature, but the principal causes 

have not yet been shown analyti cally. 

The cuspy density profile of lPrugniel fc Simienl (|l997h . 
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where x = r/r s and r s : 
profile approximately. The 
p_« 1.0 



fits the 
of the 



Sersic light 
inner cusp , 



a R c , 
index 

0.6097/n + 0.05463/n 2 (|Lima Neto et al.lll999l ; 
iMarquez et~ai1 120001 ) . Expressed in terms of incomplete 
gamma functions, the mass and moment of inertia within 
some radius are 



and 



m = 4irnb n[p 'e p s r B T[n(3 — p),bx 
/= ^ fe "M) eW r[n(5-p),bx 1/n ] ■ 
Evaluated as x — > oo, the core lever radius is finite, 
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For n = 2 we have Ri tts 3.421 r a , and for n = 4 we have 
Ri ~ 7.418 r s . The gravitational radius is finite but needs 
some numerical integration; R w /r s decreases with increasing 
n. The a 2 profile emerges from integration of a hydrostatic 
or Jeans equation, and it can be shown that the King radius 
vanishes, Rk —> as r — > 0. The dependence of the virial 
radius on (p s , r B ) is shown by grey curves in Figure |P4l 

The density index drops with radius. For all realistic 
n, i?2 exists uniquely, along with all higher slope-radii. In 
general, for index —j, we have Rj/R e — [n(j — p)/b] n . For 
n > 1 the consecutive integer-slope radii spread apart (ver- 
tical distribution of bold lines in Figure |D6|) . In contrast, 
finite polytropic halos have slope-radii spaced at shrinking 
intervals, converging at the true surface R. For n > 2, the 
rotation curve peaks inside the effective radius, R Q < R e . 
For n > 1.2 we have R < R3. The lever radius and slope-4 
radius are of similar magnitude, Ri « Ri, with Ri > Ri 
for n < 5. For n ) 1 we have Ri ^ -R3. These signature 
inequalities are potentially testable by gravitational lensing 
and kinematic studies in halo outskirts. 



Dll comparisons 

Theoretical, numerical and observational halo models are in 
principle testable by measuring enough of their signature 
radii, derived from three-dimensional and projected quanti- 
ties. Table ID1I characterises gasless polytropic halo models 
with various F values, and compares them to other models 
from the literature. 

The density slope radii are the major discriminants be- 
tween halo models. Isothermal and pseudo-isothermal mod- 
els have no R3 radius, while the Hubble, NFW and Burk- 
ert models have R3 — 00. Galaxy densities in clusters ap- 
pear to drop at least as steeply as r~ 3 in the outskirts 
(|Carlberg et al.lll997l ; lAdami et al.ll200lh . If this trend per- 
sists infinitely then it would fit NFW or Burkert interpre- 
tations. However, the detection of steeper slopes (by any 
technique) would call for more sophisticated models. All 
of the polytropic, King and Sersic models have finite R2 
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Figure D6. Signature radii (relative to R c ) for Sersic models 
with indices n. The curves are marked as in Figure |P2l In these 
halo models, the signature radii splay out enormously as n rises. 



and 7?3, set in ratios depending on F, c and n respec- 
tively. In Sersic models, the consecutive slope radii spread 
widely apart, but in the polytropes and King models these 
radii converge. Low-concentration King models have some- 
what evenly spread values of {R2, Rz, Ri, ■ ■ ■} but for high 
concentrations R2 R3 ~ R c , and so the halo could 
be mistaken for an isothermal or NFW shape if it were 
probed at intermediate radii only. For finite polytropic ha- 
los, R2/R, Rs/R, R4/R, Rw/R and R /R remain similar to 
each other in order of magnitude (even as F — > 10), while 
the projection scales Rs/R and R e /R shrink slower with 
increasing F. 

As their I and Ri/R values show, the concentrations 
and rotational properties of Hubble, NFW, Burkert and 
Hernquist models are hard to define. In one sense the mass 
is centrally concentrated. On the other hand, the outskirts 
dominate I. Therefore the ability to spin such a halo up 
or down (e.g. in tidal interactions between unbound neigh- 
bours) depends on an ad hoc truncation. The unphysical 
inability of simulacra to self-truncate may be part of the 
"angular momentum problem" of simulated galaxy forma- 
tion. The polytropic, King and Sersic models are finite and 
consistent with respect to I. 

Any real, isolated halo detaches from the Hubble flow 
possessing finite mass and energy. Its radius may also be- 
come finite, either due to intrinsic self-truncation (e.g. of a 
polytrope) or extrinsic harrassment and evaporation (e.g. 
King models). Infinite models cannot be a final or com- 
prehensive description of any real halo. They must be re- 
garded as provisional approximations only. The best phys- 
ically motivated, consistent and plausible models are the 
King, Sersic and general polytropic descriptions. It is un- 
clear how a Sersic model should adapt in a potential shared 
with other components, so it was unsuitable as the basis 
for our present study. We chose the polytropic scenario, al- 
though a generalised King model has scope to represent a 
halo suffering surface evaporation or tidal truncation. 
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Table Dl. Signature radii and masses of some representative gasless halo models. From left to right, the columns are: dark degrees 
of freedom, surface radius; inertial concentration; gravitational concentration; the density slope radii with indices —2,,— 3 and —4; the 
rotation-curve peak radius; the projected mean-light radius and effective radius; total mass. Models marked "•" and contain a 
central point mass, m* = 10~ 8 m, 10~ 6 m and 10~ 5 m respectively. 

F R R^/R. R w /R R 2 /R Rs/R Ri/R Ro/R Rs/R R e/ R m /p^l 



2 


1.253 


0.8084 


0.8000 


0.6458 


0.7286 


0.7817 


0.8733 


0.6524 


0.4636 


2.507 


3 


1.630 


0.7152 


0.7368 


0.5200 


0.6141 


0.6813 


0.7500 


0.6025 


0.3980 


3.026 


i 


2.127 


0.6222 


0.6667 


0.4126 


0.5076 


0.5830 


0.6229 


0.5500 


0.3353 


3.534 


5 


2.826 


0.5287 


0.5882 


0.3191 


0.4076 


0.4860 


0.4989 


0.4943 


0.2748 


4.040 


6 


3.891 


0.4340 


0.5000 


0.2369 


0.3135 


0.3895 


0.3810 


0.4345 


0.2162 


4.555 


7 


5.706 


0.3374 


0.4000 


0.1645 


0.2251 


0.2932 


0.2708 


0.3693 


0.1594 


5.091 


8 


9.444 


0.2376 


0.2857 


0.1010 


0.1426 


0.1965 


0.1695 


0.2955 


0.1043 


5.668 


9 


21.06 


0.1313 


0.1540 


0.04587 


0.06678 


0.09880 


0.07829 


0.2050 


0.05094 


6.323 


9.5 


44.91 


0.07286 


0.08027 


0.02164 


0.03198 


0.04947 


0.03722 


0.1444 


0.02509 


6.707 


9.5 * 


44.91 


0.07286 


0.01019 


0.02164 


0.03198 


0.04947 


0.03722 


0.1444 


0.02506 


6.707 


7 • 


5.706 


0.3374 


0.3999 


0.1645 


0.2251 


0.2932 


0.2708 


0.3693 


0.1594 


5.091 


8 • 


9.444 


0.2376 


0.04766 


0.1010 


0.1426 


0.1965 


0.1695 


0.2955 


0.1040 


5.668 


9 . 


21.06 


0.003432 


1.429("5) 


1.702(-4) 


3.674(-4) 


9.730(-4) 


2.451(-4) 


0.02658 


2.754(-4) 


6.323 


9.5 • 
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